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We model massive dark objects at centers of many galaxies as Schwarzschild black hole lenses and 
study gravitational lensing by them in detail. We show that the ratio of mass of a Schwarzschild 
lens to the differential time delay between outermost two relativistic images (both of them either 
on the primary or on the secondary image side) is extremely insensitive to changes in the angular 
source position as well as the lens-source and lens-observer distances. Therefore, this ratio can 
be used to obtain very accurate values for masses of black holes at centers of galaxies. Similarly, 
angular separations between any two relativistic images are also extremely insensitive to changes in 
the angular source position and the lens-source distance. Therefore, with the known value of mass 
of a black hole, angular separation between two relativistic images would give a very accurate result 
for the distance of the black hole. Accuracies in determination of masses and distances of black 
holes would however depend on accuracies in measurements of differential time delays and angular 
separations between images. Deflection angles of primary and secondary images as well as effective 
deflection angles of relativistic images on the secondary image side are always positive. However, 
the effective deflection angles of relativistic images on the primary image side may be positive, zero, 
or negative depending on the value of angular source position and the ratio of mass of the lens to 
its distance. We show that effective deflection angles of relativistic images play significant role in 
analyzing and understanding strong gravitational field lensing. 
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I. INTRODUCTION 

Light deflection in weak gravitational field of 
Schwarzschild spacetime is well-known since 1919 [H, 0, 
and it serves as the starting point to the learning of 
gravitational lensing (GL) theory even now [!, 0, @. 
However, light deflection in strong gravitational field of 
Schwarzschild spacetime was not studied until around 5 
decades ago Darwin [|| pioneered a theoretical research 
on GL resulting from large deflection of light in the 
vicinity of photon sphere of Schwarzschild spacetime. A 
few years later, Atkinson Q extended these studies to a 
general static spherically symmetric spacetime. In fact, 
apart from a few activities, the subject of strong gravita- 
tional field lensing remained in almost a dormant stage 
until toward the end of the last century, and this was 
due to two main reasons. First, Darwin's calculations 
showed that the images are very demagnified and there- 
fore those are very difficult to be observed with the avail- 
able observational facilities. Secondly, the known gravi- 
tational lens equation (see in [1,01) was not adequate for 
the study of lensing due to large deflection of light. As 
astronomical techniques are improving fast, it may be 
possible to overcome the observational obstacles in fu- 
ture. Therefore, an adequate lens equation was required 
for this purpose. To this end, Frittelli and Newman Q 
obtained an exact lens equation that is applicable to ar- 
bitrary spacetimes; however, this equation is difficult to 
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use in general. Further, under some physically realis- 
tic assumptions, we Q obtained a simple lens equation 
that allows arbitrary small as well as large light deflec- 
tion angles. Later, Frittelli et al. and Kling et al. [l(| 
carried out comprehensive comparative studies of the ex- 
act lens equation with our lens equation for the case of 
Schwarzschild spacetime. They found that our lens equa- 
tion works remarkably well as both approaches to gravi- 
tational lensing yield extremely close results even for light 
rays which have large deflections in strong gravitational 
field and go around the lens several times before reach- 
ing the observer. As our lens equation is easy to use and 
yields very close results to those obtained by using an ex- 
act lens equation, our lens equation has been most widely 
used in the literature for studying strong field gravita- 
tional lensing. Perlick [ll[, in a recent brief review, called 
our lens equation an almost exact lens equation. In last 
8 years, there has been a growing interest in studying 
weak as well as strong field lensing by black holes, naked 
singularities, wormholes, and some other exotic objects 

(see [n, m (n m tn [ni gj, m m us, m s m and 

references therein). 

In 0], we modeled the massive dark object (MDO) at 
the Galactic center as a Schwarzschild black hole lens, 
and, using our lens equation, studied point source GL 
due to light deflections in weak as well as strong gravita- 
tional fields. By solving the lens equation numerically, we 
obtained image positions and their magnifications. Like 
in Darwin's paper, our computations showed that, in ad- 
dition to the weak field primary (also called direct) and 
secondary images, there are theoretically an infinite se- 
quence of very demagnified images on both sides of the 
optical axis; we named them relativistic images. With 
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increase in the value of angular source position (mea- 
sured with respect to the optical axis), magnifications of 
relativistic images decrease with much faster rates com- 
pared to magnifications of primary and secondary im- 
ages. Therefore, relativistic images are not just very de- 
magnified, but are also transient. However, despite sev- 
eral other observational difficulties (discussed in and 
also in the last section of this paper), if these images 
were ever observed, we showed that it would give an up- 
per bound to the compactness of the MDO. Therefore, it 
would push black hole interpretation of the MDO at the 
Galactic center. Observation of relativistic images would 
be undoubtedly a landmark discovery in astronomy; how- 
ever, experimentalists and astronomers have to pay the 
price for these very important observations. When the 
source, the lens, and the observer are aligned, in addi- 
tion to an Einstein ring, we get theoretically an infinite 
sequence of rings; we [9j termed them relativistic Einstein 
rings. 

The central thread in this paper is a comprehensive 
study of relativistic images of Schwarzschild black hole 
lensing. The purpose of this is to theoretically investi- 
gate if possible observations of these images, compared to 
primary and secondary images, could provide more valu- 
able and accurate information about the lens. Though 
the primary and secondary images of Schwarzschild black 
hole lensing are well discussed in the literature, we also 
include these in this paper for 3 reasons. First, we do 
not take either weak or strong field approximations in 
our computations. Thus, our results are more accurate 
than known weak field limit approximate analytical ex- 
pressions could provide for lensing observables for these 
images. Therefore, these more accurate results could 
be useful for observations in near future with advanc- 
ing astronomical facilities. Secondly, it is useful to have 
thorough studies of relativistic as well as primary and 
secondary images due to a gravitational lens (with the 
same mass and distance) in the same paper so that one 
can immediately compare properties of these images. For 
observations of relativistic images and measurements as- 
sociated with them, a detailed information about their 
primary and secondary images are helpful. Thirdly, dif- 
ferential time delay of secondary image with respect to 
the primary image is though well discussed in the litera- 
ture, studies of their individual time delays were not paid 
enough attention for some reasons (discussed later in this 
paper). Therefore, we study these and obtain some im- 
portant results. 

This paper is organized as follows. In Section II, 
we first discuss a lens equation applicable to weak as 
well as strong gravitational field regions, and then give 
a brief review of deflection angles and time delays of 
light rays traveling in Schwarzschild spacetime. In Sec- 
tion III, we show that a schematic diagram for effective 
deflection angles of relativistic images give deep insight 
and provide some important information before compu- 
tations. In Section IV, we model our Galactic MDO as 
the Schwarzschild black hole lens and study variations 



in image positions, magnifications, and time delays of 
primary and secondary as well as relativistic images with 
changes in angular source position and lens to source dis- 
tance. We also study the variation of deflection angles 
for primary and secondary images and effective deflection 
angles of relativistic images with respect to changes in the 
angular source position and lens to source distance. Com- 
putations of effective deflection angles provide geometri- 
cal beauty of strong gravitational field lensing and an 
analysis of those support the reasonableness of other re- 
sults obtained through numerical computations. In Sec- 
tion V, we model MDOs at centers of many galaxies as 
Schwarzschild black hole lenses and study the variations 
in the same physical quantities (as studied in Section IV) 
with respect to changes in the ratio of mass of lens to its 
distance and the lens-source distance for a fixed value of 
angular source position. In Section VI, we discuss and 
summarize the results. 

Bozza et q/.(l7| obtained approximate analytical ex- 
pressions for image positions and magnifications of rela- 
tivistic images. Bozza and Mancini [18| further derived 
approximate analytical expressions for differential time 
delays among relativistic images. As the aim of our pa- 
per is to present accurate results and to derive some im- 
portant astrophysical implications for those, we do not 
digress to compare and contrast our results with approxi- 
mate analytical results given in [l7], EH- In Section IV, we 
briefly discuss that there are small (but significant) per- 
centage differences in results for image positions. How- 
ever, percentage differences in results for magnifications 
of images are very large. In Section V, we compare their 
approximate results with ours for differential time delays 
between two outermost relativistic images. We show that 
there are qualitative as well as large quantitative differ- 
ences. We justify our results with some arguments. 

There are some fascinating results in this paper. The 
most important among those is that relativistic images 
would provide very accurate values for masses and dis- 
tances of MDOs at centers of galaxies. The ratio of the 
mass of a Schwarzschild lens to the differential time de- 
lays between two outermost relativistic images (both of 
them either on the primary image side or on the sec- 
ondary image side from the optical axis) is almost a con- 
stant; i.e., this ratio is extremely insensitive to changes 
in the angular source position, the observer-lens distance, 
and the lens-source distance. Therefore, observation of 
relativistic images and measurements of their differential 
time delays would give very accurate values for masses 
of MDOs. Another very useful property of relativistic 
images is that variations in their angular separations due 
to changes in the angular source position and the lens- 
source distance are extremely small. Therefore, once we 
have accurate values for masses of MDOs, measurements 
of angular separations between relativistic images would 
give very accurate results for distances of those MDOs. 
Similarly, we also show that the measurement of ratio of 
fluxes of outermost relativistic images (one on each side 
of the optical axis) would help us obtain very accurate 
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result for distance of the source. 

As in our previous papers on GL [ij [l9|, [2(| [Hj], we 
use geometrized units (i.e., the gravitational constant 
G = 1 and the speed of light in vacuum c — 1, so that 
M = MG/c 2 ) throughout this paper. However, we fi- 
nally present time delays and differential time delays of 
images in the unit of minute. 



II. LENS EQUATION, DEFLECTION ANGLE, 
AND TIME DELAY 



Assuming that the angular position of source of light 
is small and the source as well as the observer are situ- 
ated at large distances from the lens (deflector), we @ 
obtained a new gravitational lens equation that allows 
for arbitrary large as well as small deflections of light. 
(The first assumption does not hurt applicability of the 
lens equation, as GL is usually meaningful only for small 
angular source positions.) The lens equation is given by 



tan (3 = tan 9 — V [tan 9 + tan (a — 9)] 



with 



V 



Dds 
D s 



(1) 



(2) 



Angular positions of an unlensed source and images are 
measured from the optical axis (the line joining the ob- 
server and the center of mass of the gravitational lens), 
and are represented by symbols (3 and 9, respectively. 
These angles when measured in clockwise and anticlock- 
wise directions from the optical axis are assigned positive 
and negative signs, respectively, a represents the total 
angle by which the light ray is deflected in the gravita- 
tional field of the lens while traveling from the source to 
the observer. Null geodesies which are bent toward and 
away from the lens have, respectively, a > and a < 0. 
Dd, Dd s , and D s stand, respectively, for observer- lens, 
lens-source, and observer-source angular diameter dis- 
tances. The values of parameter T> mathematically lie 
in the interval (0,1); however, for the lens equation to 
hold good, the value of V should not be taken too close 
to 0. The perpendicular distance from the center of mass 
of the lens to the tangent to the null geodesic at the 
source position is (see Fig. 1 in 0) 



J = Da sin 9 



(3) 



and is called impact parameter. 

The magnification of an image formed due to GL is de- 
fined as the ratio of the flux of the image to the flux of the 
unlensed source. However, according to Liouville's theo- 
rem, the surface brightness is preserved in GL. Therefore, 
the magnification p of an image formed due to gravita- 
tional lensing turns out to be the ratio of the solid angles 
of the image and of the unlensed source made at the ob- 
server [H, [j, Q . Thus, for a circularly symmetric GL, the 



magnification p of an image is obviously expressed by 

/! = mn r) (4) 

where the tangential magnification p t and the radial mag- 
nification p r are, respectively, expressed by 



sin (3 \ 1 



sin f 



and p r = 



d6 



(•5) 



Tangential critical curves (TCCs) and radial critical 
curves (RCCs) are, respectively, given by singularities in 
Pt and p r in the lens plane. However, their correspond- 
ing values in the source plane are, respectively, termed 
tangential caustic (TC) and radial caustics (RCs). The 
parity of an image is called positive if p > and negative 
if p, < 0. Sometimes terms magnifications and absolute 
magnifications of negative parity images are used synony- 
mously. If the angular source position [3 = (i.e., when 
the source, the lens, and the observer are aligned), there 
may be ring shaped image(s) [called Einstein ring(s)]; 
these images are assigned 0— parity. Note that (3 — 
does not always give Einstein ring(s) (for examples, see 

in [H^iH). 

In this paper, we thoroughly study gravitational lens- 
ing due to Schwarzschild black holes, which exterior grav- 
itational field is described by the line element 



ds 1 



2M 



df 



dr 2 



(d# 2 



sin' 



®d(j) 2 ) 



(6) 



where the real constant parameter M is the ADM 
mass. The radii of event horizon and photon sphere 
of a Schwarzschild black hole are given by R e h = 2M 
and R ps = 3M, respectively. R e h is also called the 
Schwarzschild radius. 

The bending angle a for a light ray with the closest 
distance of approach r a is given by |25|) 



a (r a ) = 2 



dr 



2M 



- 1- 



2M> 



(7) 

and the impact parameter J of the light ray is expressed 
by 



J (To) 



2M 

To 



-1/2 



(8) 



Defining a dimcnsionlcss radial distance p in terms of the 
Schwarzschild radius 2M by equation 



ft 



r 
2M 



(9) 



(for r — r , p = p ), we 0, [l9| expressed the deflection 
angle a (p) and the impact parameter J (p), respectively, 
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by 



a(p ) 



and 



dp 



-1/2 



J(p )=2Mp 1 



(10) 



(11) 



For computations of magnifications of images, the first 
derivative of deflection angle a with respect to is 
needed, which is given by H,[H[ 



da 



(12) 



where the first and second factors on right side of this 
equation are, respectively, given by 



a' (p ) = 



2p 



(4p-3)dp 



(3-2p) a ^(-) a (l-^)-(l-j) 



(13) 



r 



and 



dp. 



It:(2Po-3) 



(14) 

Time delays for images of gravitational lensing are 
given by 3 terms: the first and second terms with pos- 
itive sign are, respectively, the travel time of the light 
from the source to the point of closest approach and 
from that point to the observer, and the third term 
with a minus sign is the light travel time from the 
source to the observer in the absence of any gravita- 
tional field. Solving null geodesic equations for general 
static spherically symmetric spacetime, Weinberg in his 
classic book (25| obtained the time required for light to 
travel from a source at coordinates {r, 1?, tt/2, if = <fi} 
to the closest point of approach (to the lens) at coordi- 
nates {ro, 7r/2, ip = ifz}- Using this result, time delay 
of images of Schwarzschild lensing can be expressed as 
(see Eqs. (23)- (25) in [U) 




A, 



dp 



AV, 



° dp 

f(P) ' J p f(P) 



-D s sec (3 (15) 



j±) tan* ). 



Xo = —, (16) 
2M' v ' 



and 



f(p) 




1 

1 - - 
P 



1 

Po 



Time delay of a gravitationally lensed image may be in 
general positive, zero, or negative; for examples, see in 



our paper [2l| . However, time delays are always positive 
for images of Schwarzschild lensing. 

It is worth mentioning that Eq. (4.67) in a classic 
book on GL by Schneider et al. [3| gives time delays of 
gravitationally lensed images. That equation contains an 
additive constant term. The authors clarified that the 
constant term is the same for all rays from the source 
plane to the observer. Therefore, this term cancels for 
computations of differential time delay between 2 im- 
ages. However, as the value for the constant term is not 
obtained, that equation cannot be used to compute time 
delays of images. This is why we follow the approach 
given in Weinberg's book (2f| and we discussed that in 
this section. 



III. EFFECTIVE DEFLECTION ANGLES OF 
RELATIVISTIC IMAGES 

It is important to first discuss in brief a few new terms 
we defined in our previous paper [9( . Then, we will show 
that these definitions with some arguments, reveal ge- 
ometrical beauty of strong field Schwarzschild lensing. 
This also helps predict some results without computa- 
tions and thus provides consistency check for results ob- 
tained through numerical computations. 

If a lens is very compact, then a light ray passing close 
to it will suffer a large deflection and therefore will loop 
around the lens once, twice, thrice, or many times (de- 
pending on the closest distance of approach from the 
center of the lens) before reaching the observer. We 
defined relativistic images of GL as those images which 
occur due to light deflections by angles a > 3ir/2. Sim- 
ilarly, for the angular source position /3 — 0, we defined 
relativistic Einstein rings as those ringed-shaped images 
which can form due to light deflections by angles a > 2ir. 
Relativistic Einstein rings are thus relativistic images for 
the case of (3 = 0. It is useful to define order of relativis- 
tic images on each side of the optical axis. We assign 
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the order 1 to the outermost relativistic images on both 
sides of the optical axis, 2 for adjacent inner ones, and 
so on. Thus, according to this definition, the outermost 
relativistic Einstein ring also has the order 1. 

In a recent paper [2l| . we discussed that the existence 
of a photon sphere enclosing a lens is a sufficient (not nec- 
essary) condition for the formation of relativistic images. 
A sufficiently compact lens can give rise to relativistic 
images even if the lens is not covered inside a photon 
sphere. Therefore, the lens need not be a black hole to 
produce these images. 

In Q , we defined a term effective deflection angle of a 
relativistic image, which we now express as follows: 

a e (po) = a (p ) - 2mr, (18) 

where n is a positive integer that represents the number 
of loops (turns) a light ray makes around the lens before 
reaching the observer, and the superscript e on a stands 
for the word effective. (In fact, the above equation can be 
also applied to primary and secondary images, because 
n = correspond to those images. However, we prefer 
to call those as known in the literature instead of calling 
them images of 0-order.) a(po), as given by Eq. (TIT)]) , 
is the usual (total) deflection angle for the light ray with 
the closest scaled distance of approach po making n loops 
around the lens. Thus, similar to the cases of primary 
and secondary images, the effective deflection angle of a 
relativistic image turns out to be the angle between the 
tangents on the null geodesic at the source and at the 
observer positions. In the following, we now introduce 
new symbols we use for physical quantities associated 
with primary, secondary, and relativistic images. 

New symbols. — We use subscripts p and s, respectively, 
for primary and secondary images. Similarly, subscripts 
np and ns (n — 1, 2, 3, ...) stand, respectively, for the im- 
ages of order n on the same side as the primary and sec- 
ondary images. For example, a\ p and af 8 , respectively, 
stand for effective deflection angles of relativistic images 
of order 1 on the same side as the primary and secondary 
images. The same applies to symbols for angular image 
positions, magnifications, and time delays. 

It is well-known that, on the same side as the source 
from the optical axis, the Schwarzschild black hole lens- 
ing gives rise to the primary image which is formed due to 
light deflection in weak gravitational field without loop- 
ing of the ray of light around the lens. On the other 
hand, relativistic images on the same side as the source 
are produced due to looping of light rays around the lens, 
which is caused by large deflection angles a > 3tt/2 in 
strong gravitational field. It is natural to ask if, exclud- 
ing the primary image, there is any other image (on the 
same side as the source) which can form without looping 
of the light ray around the black hole. Computations give 
no such solutions to the lens equation. Therefore, on the 
same side as the source, there is only one (i.e., the pri- 
mary) image which forms due to light deflection in weak 
field without looping of the light ray around the lens and 
there are relativistic images which arise due to looping 



of light rays around the lens in strong gravitational field. 
In the following, we will show that a simple geometrical 
argument beautifully supports these numerical results. 

See Fig. 1. Assume that, on the same side as the 
source from the optical axis, two light rays emitted from 
the source S reach the observer O without looping around 
the lens. SP± and SPz are tangents on 2 null geodesies 
at the source position, and, similarly, P\0 and P2O are, 
respectively, tangents on those null geodesies at the ob- 
server position. &i and A2 are deflections angles corre- 
sponding to 2 light rays we considered and 71 and 72 are 
their respective supplementary angles. The schematic di- 
agram (right of Fig. 1) shows that 72 > 71. This implies 
that a.2 < &i, which is not allowed according to the a 
vs J/M graph (see the left of Fig. 1); i.e., a decrease 
in the impact parameter should increase a. Thus if the 
light path SP\0 is allowed, then SP2O is not allowed 
and therefore we conclude that there can be only the pri- 
mary image on the same side as the source without the 
light ray looping around the lens. A similar argument 
also demonstrates that there can be only one image (i.e., 
secondary) on the opposite side from the source without 
a light ray going around the lens. 

Now consider that SP2 and P2O represent tangents, 
respectively, at the source and the observer position on 
a null geodesic that loops around the lens once before 
reaching the observer. Therefore, 612 — c\\ p < &i, which 
also reflects in numerical computations in next section. 
For a given value of /3, the schematic diagram also shows 
that the effective deflection angle of relativistic images 
on the same side as the primary image decrease with in- 
crease in its order. Similarly, for a given angular source 
position, the effective deflection angles of relativistic im- 
ages (on the same side as the secondary) decreases with 
the increase in its order. These conclusions based on 
simple geometrical analysis are reflected in results of our 
numerical computations in the next section. 

In Fig. 2, we show that effective deflection angles of 
relativistic images on the same side as the primary image 
can be positive, zero, or negative depending on the value 
of the angular source position. However, for relativistic 
images on the same side as the secondary image, effec- 
tive deflection angles are always positive. Consider the 
first order relativistic images on both side of the optical 
axis; i.e., one on the primary image side and the other 
on the secondary image side. SC\ and C\0 are, respec- 
tively, tangents on null geodesies (giving rise to the 1 st 
order relativistic image on the primary image side) at 
the source and observer positions; C± is their point of 
intersection. Similarly, SC2 and C2O are, respectively, 
tangents on null geodesies (giving rise to the I s * order 
relativistic image on the secondary image side) at the 
source and observer positions; C2 is their point of inter- 
section. For small angular source position 0, the effective 
deflection angle a\ > (see the extreme left diagram). 
As the value of (3 increases, the value of af decreases 
to zero value (see the second and third diagrams from 
the left). When df = 0, the angular source position 
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Source Relativistic image Primary image 




FIG. 1: (color online). Left: The deflection angle a is plotted against the dimensionless scaled impact parameter J/M. The 
arrow attached to the curve indicates that a decreases with increase in the value of J/M. Right: S, L, and O represent 
positions of the source, the lens, and the observer, respectively. SPi and SP2 are tangents on 2 null geodesies at the source 
position, whereas P\0 and P2O are tangents on, respectively, the same pair of null geodesies at the observer position. d?i 
and o?2 are light bending angles, whereas 71 and 72 are their respective supplementary angles. /3 and 8 P stand, respectively, 
for the angular positions of the source and the primary image. Angles in this schematic diagram are greatly exaggerated. 
72 > 71 =>■ &2 < &\. According to the a vs J/M plot (on left side), 612 < &\ is possible only if SP2 and P2O are, respectively, 
tangents at the source and the observer positions on a null geodesic which loops around the lens at least once giving rise to a 
relativistic image. 



of this relativistic image and the source coincide. We 
denote this critical value of the angular source position 
as I3\ c . (subscript c stands for the word critical and 1 
stands for the 1 st order relativistic image.) A further 
increase in the value of (3 makes the effective deflection 
angle 6l\ p < and its value keep decreasing with increase 
in the value of (3. On the other hand, the effective de- 
flection angles of relativistic images on the same side as 
the secondary image increases with the increase in the 
value of (3. These important conclusions based on simple 
geometrical analysis also appear in our numerical results 
in the next section. Therefore, the geometrical analysis 
using Fig. 1 and Fig. 2 also supports correctness of our 
computations. 



IV. GRAVITATIONAL LENSING BY THE 
GALACTIC MDO 

In this section, we model the MDO at the center of 
our Galaxy (the Milky Way) as a Schwarzschild black 
hole lens and study point source GL in a great detail. 
The MDO has the mass M = 3.61 x 1O 6 M and is at 
the distance D c i = 7.62 kpc from us [2(|. Therefore, 



M/D d « 2.26 x 10" 11 (note that M = MG/c 2 ). In a 
recent paper [2l| . we already obtained angular positions, 
deflection angles, magnifications, and time delays for pri- 
mary and secondary images for several values of angular 
source position (3 for T> = D c i s /D s — 0.5 (i.e., when the 
lens is situated halfway between the observer and the 
source). We also computed differential time delays of 
secondary images with respect to primary images. For 
comparison and continuity in discussion, we use those 
results in this paper and also put those in Table I in 
the Appendix of this paper. Now considering D = 0.5 
and using Mathematica, we numerically solve the grav- 
itational lens equation (TT]) for a large number of values 
for (3 and obtain image positions for first and second or- 
der relativistic images on both sides of the optical axis. 
We further obtain deflection angles, magnifications, time 
delays, and differential time delays for these images. Us- 
ing Eq. ([15]). we obtain effective deflection angles for 
these relativistic images. We put these results in Tables 
II and III. Though, we computed for a large number of 
values for (3, we put only a few data in tables; however, 
we use all those for figures. We further repeat the entire 
computations for primary and secondary as well as rela- 
tivistic images for V — 0.05 and 0.005 to see the effects 
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FIG. 2: (color online). Schematic diagram showing the variation in effective deflection angles of relativistic images with respect 
to increase in the value of angular source position /3. S, L, and O stand for the position of the source, the lens, and the 
observer, respectively. SCi and SC2 are tangents on 2 null geodesies at S, and CiO and C2O are tangents on corresponding 
null geodesies at O. a\ p and df s stand for effective deflection angles of relativistic images of order 1, respectively, on the 
primary and secondary image sides; 9i p and 9i s are their respective angular image positions. Angles are greatly exaggerated. 



of changes in image positions, deflection/effective deflec- 
tion angles, magnifications, time delays, and differential 
time delays due to change in the lens-source distance. 
(With the observer-lens distance Dd fixed, a decrease 
in the value of V decreases the source-lens distance.) 
Throughout our computations in this paper, we never 
take either weak or strong gravitational field approxima- 
tion and therefore our computations and hence results 
are exact in this sense. In the following paragraphs, we 
will now discuss results for GL by the Galactic MDO. 

In Fig. 3, we plot the (absolute) angular positions 
(measured from the optical axis) of primary and sec- 
ondary images, and their separations against the angular 
source position (3 for T> = 0.5,0.05 and 0.005. As it 
is well-known that, for a given value of V, the angular 
positions of primary and secondary images, respectively, 
increase and decrease with an increase in the value of (3. 
For a given value of /?, the angular positions of primary as 
well as secondary images increase with an increase in the 



value of T>. The angular radius of Einstein ring increases 
with increase in the value of T>. The angular separation 
between primary and secondary images increases with 
increase in f3 and T>. The angular source positions of 
relativistic images are very insensitive to changes in the 
values of (3 and T>. The angular position of the first order 
relativistic images on both sides of the optical axis have 
extremely close values; however, 9\ p > \9i s \ for all values 
of (3, excluding, of course, at (3 — for which 9\ p = \6\ s \. 
The same is true for any pair of second or higher order 
relativistic images. As 9 np and \9 ns \ (for the same value 
of n) have extremely close values, we plot image posi- 
tions of relativistic images only on the same side as the 
primary image. For theoretical interest, it is worth in- 
vestigating variation in the value of 9 np and \9 ns \ with 
changes in the value of (3 and T>, though variations are 
extremely small. For relativistic images of a given order 
n and for a given value of V, the values in 9 np and \9 ns \, 
respectively, increase and decrease with the increase in 
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FIG. 3: (color online). Top left and middle: The angular positions of primary images 9 P , secondary images \6 S \, and their 
separations 9 P — 9 S are plotted against the angular source position /3 for T> — 0.5,0.05 and 0.005. Top right: The angular 
positions of relativistic images (on the same side as the primary image) of the first order 9\ p and of the second order 02 P 
are plotted against /3 for the same values of T> as in the figures on left. The curves for 8i p (for different values of T>) 
intersect for /3i c ~ 24.3028 /ias, whereas those for 02 P intersect for fiic ~ 24.2724 \ias. Below: The angular separations among 
relativistic images vs the angular source position are plotted for T> = 0.5,0.05 and 0.005. 9 np and 9 na (n — 1,2) stand 
for angular positions of relativistic images on the primary and the secondary image sides, respectively. The Galactic MDO is 
modeled as the Schwarzschild lens, which has mass M = 3.61 x 1O 6 M0 and is situated at the distance Dd = 7.62 kpc so that 
M/D d « 2.26 x 10~ n . 



the value of (3, though the increase and decrease are ex- 
tremely small. However, their dependence on T> is much 
more fascinating. For (3 — or a small value, 9\ p is 
greater for a greater value of T>. As j3 increases, there is 
situation when Q\ v is the same for all values of V. At this 
critical angular source position f3\ c ss 24.3028 microarc- 
sec (/ias), the effective deflection angle af p — for all 
values of T>. For a further increase in the value of /3, 9\ p 
is smaller for a greater value of T>. The same happens for 
any relativistic images on the same side as the primary 
image; however, critical angular source position [3 nc de- 
creases with an increase in the order n of the image; for 
example, foe ~ 24.2724 /ias. On the other hand, for any 
fixed value of (3, image positions of relativistic images 



(on the secondary image side) \0 ns \ always increase with 
an increase in T). The angular separation between rela- 
tivistic images of the first order (i.e, 8i p — 0i s ) increases 
with an increase in the value of (3 (for a fixed V). The 
increase rate of the angular separation with increase in (3 
is higher for lower value of V. For (3 = or a small value, 
this angular separation is higher for a higher value of T>. 
(This is qualitatively similar to the case of image separa- 
tion between primary and secondary images.) However, 
for large value of (3, the angular separation # lp — ls is 
higher for lower value of T>. For any two values of T>, 
there exist a (3 at which values of 9\ p — 9i s are equal for 
both T>. The angular separations between outermost 2 
relativistic images (both of them either on the primary 
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FIG. 4: (color online). Top left: The deflection angles of primary images a p and secondary images & s are plotted against the 
angular source position j3 for T> = 0.5, 0.05 and 0.005. Top middle and right: The effective deflection angles of the first order 
relativistic images on the same side as the primary image a\ v and on the same side as the secondary image a\ s are plotted 
against /3 for T> = 0.5, 0.05 and 0.005. Below left and right: The effective deflection angles of the second order <k% v (left) and 
of the first order a\ v (right), both on the same side as the primary image, are plotted against /3 for T> — 0.5, 0.05 and 0.005 in 
the vicinity of zero effective deflection angle. /?2c ~ 24.2724 fias and /3i c ~ 24.3028 jj,as, where {3\ c and (3\ c are, respectively, 
critical angular source positions for the second and first order relativistic images. The gravitational lens is the same as for the 
Fig. 3. 



or on the secondary image side) have just opposite qual- 
itative dependence on (3 and T> in the following sense. 
dip — &2p increases with increase in (3 (for a fixed £>), 
but decreases with increase in V (for a fixed (3). On 
the other hand, \9\ s — #2s decreases with increase in f3 
(for a fixed T>), but increases with increase in T> (for a 
fixed (3). The increase/decrease rate with change in the 
value of (3 (for a fixed T>) is smaller for higher value of 
T>. The variations in angular separations between rela- 
tivistic images are very small with respect to changes in 
(3 and T>. Among angular separations 9\ p — 9i s , 9\ p — 02 P , 



and \9i s — #2 S |, the first one is the least sensitive to those 
changes. As relativistic images would be observationally 
important for very small values of (3, we conclude that 
angular separations between relativistic images of our in- 
terest are extremely insensitive to change in the value of 
V. 

In Fig. 4, we show variations in deflection angles for 
primary and secondary as well as effective deflection an- 
gles for relativistic images with respect to changes in val- 
ues for (3 and V. It is known that for a given value of 
2?, the deflection angles for primary and secondary im- 
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FIG. 5: (color online). Top left and middle: The magnifications of primary images fi p , the absolute magnifications of secondary 
images \fi s \, and their ratios \fi p /fj, s \ are plotted against the angular source position (3 for different values of D. Top right and 
below left : The magnifications of relativistic images (on the same side as the primary image) of the first order pi p and the 
second order \xi v are plotted against the angular source position /3 for D = 0.5,0.05 and 0.005. Below middle and right: The 
magnifications ratios |/ii p //ii s | (where fii s stands for magnification of the first order relativistic image on the same side as the 
secondary image) and fJ,i p /fi2 P vs {3 are plotted for the same values of T> as in the figure on below left. The lens is the same as 
for the Fig. 3. 



ages, respectively, decrease and increase with an increase 
in the value of (5. For a fixed value of /?, deflection angles 
for these images increase with decrease in the value of V. 
Like primary and secondary images, the deflection angles 
of relativistic images are always positive, and the same 
is true for effective deflection angles of relativistic images 
on the secondary image side. However, for a fixed value 
of M/Dd, the effective deflection angles of relativistic im- 
ages on the same side as the primary image decrease with 
an increase in the value of /3, and can be positive, zero, 
or negative depending on the value of /?. For a relativis- 



tic image (of any order) on the same side as the primary 
image, there exists a critical value of f3 for which the de- 
flection angle is zero. These results are as expected from 
the schematic diagram in Fig. 2. Our numerical compu- 
tations give (3 2c « 24.2724 fias and /3 Xc « 24.3028 jias 
showing that /32c < Pic- The effective deflection angle is 
positive, zero, and negative, respectively, for the angular 
source position less than, equal to, and greater than the 
critical value of the angular source position. The critical 
source positions are independent of the value of T>. As 
expected from the schematic diagram (see Fig. 2), our 
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FIG. 6: (color online). Top: The time delays of primary images r p and secondary images r s , and the differential time delays 
of the secondary images with respect to their respective primary images (i.e., t s — t p ) are plotted against the angular source 
position j3 for T> = 0.5, 0.05 and 0.005. Below: The time delays of relativistic images (on the same side as the primary image) 
of the first order rip, second order T2 P , and the differential time delays (i.e., T2 P — Tip) axe plotted against j3 for the same values 
D as in figures on top. The lens is the same as for the Figs. 3 through 5. 



computations also show that effective deflection angles 
of relativistic images on the secondary image side are al- 
ways positive and increase with increase in the value of 
[3 (for fixed V). 

In Fig. 5, we show changes in (absolute) magnifications 
of primary and secondary images as well as relativistic 
images with changes in the values of (3 and T>. Images 
on the same side as the source and opposite side from 
the source have, respectively, positive and negative mag- 
nifications and therefore have, respectively, positive and 
negative parities. The (absolute) magnifications of pri- 
mary, secondary, and relativistic images decrease with an 
increase in the value of angular source position [3. How- 



ever, there are 2 important differences between absolute 
magnifications of primary-secondary pair and relativistic 
images. First, the absolute magnifications of relativis- 
tic images are extremely small and decrease much faster 
than those of primary and secondary images with an in- 
crease in the angular source position (3. Secondly, as 
opposed to the case of primary and secondary images, 
for a fixed value of f3, the absolute magnifications of rel- 
ativistic images increase with decrease in the value of T>. 
Therefore, it would be easier to observe relativistic im- 
ages of those sources which are relatively nearer to the 
lens. The absolute magnifications of relativistic images 
of the same order on each side of the optical axis have 
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FIG. 7: (color online). The ratios of mass M of the lens to differential time delays of images are plotted against the angular 
source position /3 for T> — 0.5, 0.05 and 0.005. t p and t s stand for time delays of primary and secondary images, respectively. 
rip and T2 P , respectively, represent time delays of relativistic images (on the same side as the primary image) of orders 1 and 
2. The gravitational lens is the same as for the Figs. 3 through 6. 



extremely close values (images on the same side as the 
primary image have though slightly higher value than im- 
ages on the same side as the secondary image). This is 
why we plot only for relativistic images on the primary 
image side. For a fixed value of V, the ratios of absolute 
magnifications of primary and secondary images |/x p //x s |, 
and of relativistic images and £ti p //Z2 P increase 

with increase in the value of (3. However, for a given 
value (3, these decrease with an increase in V. Compared 
to the (absolute) magnifications ratio of primary and sec- 
ondary images, the ratios of (absolute) magnifications of 
relativistic images are much less sensitive to changes in 
j3 and V . 

In Fig. 6, we first plot time delays of primary and sec- 
ondary images, and differential time delays of secondary 
images with respect to their respective primary images 
against the angular source position (3 for V — 0.5, 0.05, 
and 0.005. For Einstein rings (/3 = case), time delay 
decreases with a decrease in the value of T>. For any given 



value of V, the time delays of primary and secondary im- 
ages, respectively, decrease and increase with increase in 
the value of (3. Similarly, for any given value of /?, time 
delay of a primary image decreases with a decrease in the 
value of T>. However, there is no such simple dependence 
for the time delays of secondary images on T>. For a small 
value of (3, time delay of secondary image is smaller for 
smaller T> and the difference decreases as (3 increases. For 
a certain value of (3, time delays for secondary image for 
2 different values of V become equal. For a further in- 
crease in f3, time delays for secondary images are higher 
for lower value of T> and the difference keeps increasing 
with increase in f3. For a fixed value of T>, the differential 
time delay of secondary image with respect to the pri- 
mary image increases with increase in (3. However, for 
a fixed value of f3, this differential time delay increases 
with decrease in T>. We now plot time delays of rela- 
tivistic images of the first and second orders (both on 
the same side as the primary image) , and the differential 
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time delay of the first with respect to the second against 
ft for V = 0.5, 0.05, and 0.005. The time delays of rel- 
ativistic images of the same order on each side of the 
optical axis have extremely close values (images on the 
same side as the primary image though have lower values 
than images on the same side as the secondary image). 
This is why we plot only for relativistic images on the pri- 
mary image side. The differential time delay — Tip) 
has simple dependence on ft and V. For a fixed value of 
T>, the differential time delay increases with an increase 
in the value of ft; however, for any fixed value of ft, the 
differential time delay increases with a decrease in the 
value of T>. For fixed ft and T>, time delays of relativistic 
images increase with increase in the order; i.e., inner rel- 
ativistic images have higher time delays relative to outer 
relativistic images. For a given value of T> and order n, 
time delay of a relativistic image increases with increase 
in the value of ft; the rate of increase is higher for lower 
values of T>. However, for a fixed value of ft and order 
n of relativistic image, the dependence of time delays on 
T> is not so simple. For ft = or a small value, the time 
delay of a relativistic image of a given order is smaller for 
a smaller value T>; however, for a certain value of ft, both 
equal and as ft increases further, time delay is higher for 
a lower value of V. 

In Fig. 7, we plot ratios of the mass of the lens to dif- 
ferential time delays among images against the angular 
source position ft for T> = 0.5, 0.05, and 0.005. We choose 
4 differential time delays: (a) differential time delay of 
secondary image with respect to the primary image; i.e., 
r s — t p , (b) differential time delay of the first order rel- 
ativistic image (on the same side as the primary image) 
with respect to the primary image, i.e., t\ p — t p , (c) differ- 
ential time delay of the first order relativistic image (on 
the same side as the primary image) with respect to the 
secondary image, i.e., t\ v — t s , and (d) differential time 
delay of the second order relativistic image with respect 
to the first order relativistic image (both on the same 
side as the primary image), i.e., r 2p — t\ p . We do not 
consider some other combinations for differential time 
delays for obvious reasons; for example, we do not use 
Tig — Ti p , because these are too small (see Tables 2 and 
3) to be measured possibly in several decades to come. 
For a fixed value of V, ratios M/(t s — t p ),M / (ti p — r p ), 
and the ratio M/(t2 P — t\ p ) decrease and M/(t\ p — t s ) 
increase with increase in the value of ft. For any given 
value of ft, M/(t s — t p ) and M/(r 2p — t\ p ) decrease and 
M/ (ti p — t s ) increases with a decrease in the value of V; 
however, dependence of M/(t± p — t p ) on T> is somewhat 
complex and fascinating. For ft — or a small value, 
the ratio Mj(j\ v — r p ) is higher for smaller V. As ft in- 
creases, 2 curves for 2 different values of T> intersect and 
hence this ratio is the same for both values of T>. For a 
further increase in ft, the ratio is now higher for higher 
value of T>. Fig. 7 shows that the ratio M/(t2 P — Ti p ) is 
the most insensitive to changes in values of ft and T>. In 
fact, as the relativistic images can be observed only for 
ft = or a very small value, the variation in the ratio 



M/(t2 P — Ti p ) due to change in T> is extremely small. In 
the next section, we will show that this ratio is in fact 
extremely insensitive to change in the value of M/Dd as 
well. Therefore, the physical quantity M / (r 2p — T\ p ) can 
be approximately considered as a constant, which can be 
used to compute very accurate values for masses of black 
holes once differential time delays r 2p — T\ p are measured. 

We mentioned in the first section of this paper that 
Bozza et al. (l7l | analytically obtained approximate ex- 
pressions for image positions and magnifications of rel- 
ativistic images. In order to calculate angular positions 
of these images, they first obtained an expression for ef- 
fective deflection angles (though they did not call that 
by this name). Here, we briefly compare our results with 
their by giving some examples. For the purpose of com- 
parison, we consider the MDO at the center of the Milky 
way as the lens. This lens has M/Dd ~ 2.26 x 10 -11 . 
We consider the lens to be situated halfway between the 
source and the observer (i.e., T> = 0.5) and the angu- 
lar source position ft = 1/ias. Compared to our results, 
Bozza et al. expressions give f» 0.5% higher values for 
each of the following: a\ (effective deflection angle of 
the relativistic image of order 1 on the primary image 
side), 6\p (angular position of the relativistic image of or- 
der 1 on the primary image side), and Q\ p — Q\ s (angular 
separation between relativistic images of the first order). 
Though, percentage differences appear small, these are 
significant for 2 reasons. We have shown that the angular 
positions of relativistic images and their separations are 
extremely insensitive to changes in the angular source po- 
sition as well as lens-source distance. In view of this fact, 
the above percentage differences are significant. In the 
next section, we show that angular separation between 
2 relativistic images can be used to obtain very accurate 
value for distance of a lens. Therefore, percentage errors 
in Bozza et al. results will decrease accuracies in deter- 
mination of distance of MDOs. Secondly, observation of 
relativistic images would provide a method to test the 
general theory of relativity against alternative theories 
of gravity in strong gravitational field region. Angular 
positions of relativistic images of the same order and the 
same parity in different theories of gravity are expected 
to be very close. Therefore, very accurate theoretical re- 
sults for image positions would be required to compare 
different theories of gravity. We now compare results for 
magnifications of relativistic images due to the same lens. 
For V = 0.5 and ft — 1/ias, Bozza et al. result yields 
w 372% higher value than our for ^ lp (magnification of 
the first order relativistic image on the primary image 
side). This very large percentage difference appears to 
be due to unrealistic drastic approximation they took in 
their calculation. Moreover, according to their result, the 
absolute magnification of relativistic images of the same 
order are equal; i.e. \^ n p/^ns \ = I , which is obviously not 
correct due to the asymmetry (ft ^ 0). Our results show 
that {fJ-np/ fJ-nsl > 1, as expected. For qualitative similari- 
ties between Bozza et al. and our results, see [l7|. Bozza 
and Mancini [l8| also analytically obtained approximate 
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expressions for differential time delays among relativistic 
images. In the next section, we show that there are again 
large percentage errors in their results. 



V. GRAVITATIONAL LENSING BY MDOS AT 
CENTERS OF MANY GALAXIES 

In this section, we model MDOs at centers of 40 galax- 
ies as Schwarzschild black hole lenses and, like in the pre- 
vious section, study point source GL by them. Gebhardt 
[27l ] tabulated updated values of masses and distances of 
many MDOs. In Table 4, we consider most of those and 
arrange in the decreasing order of M/Dd (i.e., the ratio 
of mass M and the distance Dd) of MDOs. (Only for the 
Galactic MDO, we use the updated values of mass and 
distance given in [26| .) The aim of this section is to study 
variations in angular positions of images and their sep- 
arations, deflection angles (effective deflection angles for 
relativistic images), magnifications and their ratios, time 
delays, differential time delays, and ratios of the mass of 
the lens to differential time delays due to changes in the 
value of M/D d for V = D ds /D s = 0.5, 0.05, and 0.005. 
We mentioned in the first section that the central thread 
in this paper is the study of relativistic images and we 
know that these images may be observed only when the 
lens components (the observer, the lens, and the source) 
are perfectly or highly aligned. In view of this we take the 
angular source position /3 = lfias for computations. (As 
we also want to compute magnifications of images due to 
point source GL, we do not take (3 = 0.) For (3 = 1/zas, 
V = 0.5, and different values of M/Dd for several MDOs, 
we numerically solve the lens equation to obtain image 
positions of primary and secondary as well as relativistic 
images of orders 1 and 2. Further, we compute deflec- 
tion angles, magnifications, time delays, and differential 
time delays for primary and secondary as well as rela- 
tivistic images. From deflection angles of relativistic im- 
ages, we compute effective deflection angles. In Tables 5 
and 6, we present results, respectively, for primary and 
secondary images, and for relativistic images. Though 
we do not display results for deflection angles of primary 
and secondary images, and effective deflection angles for 
relativistic images in tables, we use those in figures. Fur- 
ther, we repeat the entire computations for T> — 0.05 and 
0.005. As in the previous section, we do not take either 
weak or strong gravitational field approximation in any 
part of our computations and therefore our results are 
very accurate. With all results available, we present sev- 
eral plots and discuss these in the following paragraphs. 
We do not present some results for relativistic images on 
the secondary image side, because computations show 
that those results for the same order relativistic images 
on both sides of the optical axis are extremely close and 
therefore graphs for those do not appear resolved on fig- 
ure. 

In Fig. 8, we first plot the angular positions of primary 
and secondary images, and their separations against the 



ratios of the mass of the lens to its distance (i.e., M/Dd) 
for T> = 0.5, 0.05, and 0.005. As expected from well- 
known analytical expressions for primary and secondary 
image positions, the angular positions of these images, 
for a given value of the angular source position (3, in- 
crease with increase in the values of M/Dd and V. As 
we have taken /3 = 1/ias (a very small value), the curves 
for primary and secondary images are too close to ap- 
pear resolved on the figure. However, angular position 
of a primary image is always greater than that of a 
secondary image. The angular separation between pri- 
mary and secondary images increases with increase in 
the value of M/Dd and T>. Further, we plot position of 
the first order relativistic image 6± p against M/Dd for 
V = 0.5, 0.05, and 0.005. (We do not plot image posi- 
tions for the second order relativistic images Oip , because 
on chosen scales these do not appear separate from the 
curve for 6\ p .) For a fixed value of T> and order n, 9 np 
increases with increase in M/Dd- The curves for differ- 
ent values of V are too close to appear separate on the 
figure. Note that there is a critical value of M/Dd at 
which curves for different values of T> intersect. For the 
first order relativistic images on the primary image side, 
(M/D d )ic « 9.31854 x 10~ 13 . For M/D d < (M/D d ) lc 
and M/Dd > (M/ Dd)i c , the value of 6\ p are, respec- 
tively, lower and higher for higher value of T>. Obvi- 
ously, for M/Dd = (M/Dd)ic, 0i P is the same for any 
value of V. Similarly, the critical value for the second 
order relativistic images, (M/D d ) 2c « 9.33022 x 10~ 13 
and the above results apply for 02 P also. We find that 
(M/Dahc > (M/D d )ic- For any given value of M/D d 
and order n, the angular positions of relativistic images 
(on the secondary image side) are higher for higher T>. 
In the next paragraph, we show that at critical values 
for M/Dd, the effective deflection angles of relativistic 
images on the primary side are zero. The critical value 
of M/Dd for any give order of relativistic image depends 
on the value of the angular source position. We further 
plot angular separations between first order relativistic 
images (i.e., 9\ p — 9i s ) and between relativistic images 
(both on the primary image side) of the first and second 
orders (i.e., 6\ p — #2 P ) against M/Dd for V = 0.5,0.05, 
and 0.005. Our results show that variations in angular 
positions of relativistic images and their separations are 
extremely small (insignificant) due to change in the value 
of V. 

In Fig. 9, we first plot deflection angles for primary and 
secondary images against M/ Dd (the ratio of mass to dis- 
tance of lens) for T> = 0.5,0.05, and 0.005. As we have 
taken the angular source position [3 a very small value, 
the curves for primary and secondary images do not ap- 
pear separate on the figure. The deflection angles for 
these images increase with increase in the value of M / Dd 
(for a fixed value of T>) and decrease with an increase in T> 
(for a fixed value of M/ Dd) ■ We now plot the effective de- 
flection angles of the first order relativistic images on each 
side of the optical axis. The dependence of these effective 
deflection angles on T> and M/Dd is qualitatively similar 
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FIG. 8: (color online). Top left and middle: The angular positions of primary images 9 P (represented by continuous curves), 
secondary images \9 S \ (represented by dotted curves), and their separations 9 P — 9 S are plotted against M/Dd of MDOs at 
centers of many galaxies for T> = 0.5, 0.05 and 0.005. Top right: The angular positions of the relativistic images (on the same 
side as the primary image) of the first order 9i p are plotted against M/Dd for the same values of T) as in figures on left. 
The curves for 9\ p for different values of T> intersect for (M/Dd)ic ~ 9.31854 x 10~ 13 . Below middle and right: The angular 
separations 9\ p — 02 P and 9i p — 9\ s among relativistic images vs M/Dd are plotted for the same values of T> as in figures on 
top. #2 P and 9i 3 stand for angular positions of relativistic images of second order on the primary image side and of first order 
on the secondary image side, respectively. The angular source position /3 = lfias for all figures. 



as for deflection angles for primary and secondary images. 
However, there is an important difference. The deflection 
angles of primary ans secondary images, and effective 
deflection angles of relativistic images on the secondary 
side are always positive. However, the effective deflection 
angles of relativistic images on the primary image side 
are negative, zero, or positive depending on the value of 
M/Dd- The critical value of the ratio (where the effective 
deflection angles are zero) for the first and second order 
relativistic images are (M/D d ) lc ps 9.31854 x 10~ 13 and 
{M/D d ) 2c ~ 9.33022 x 1CT 13 , respectively. These values 
are exactly the same as we obtained for intersections of 
curves for image positions for different values of V (see 
the previous paragraph). 

In Fig. 10, we first plot (absolute) magnifications of 
primary and secondary images, and their ratios against 
M/D d for V = 0.5, 0.05, and 0.005. The magnifications 



increase with increase in the value of M/Dd (for a fixed 
valued of T>) as well as T> (for a fixed value of M/Dd). As 
the chosen angular source position is very small, the ratio 
of these magnifications is very close to 1. We then plot 
magnifications of relativistic images of first and second 
orders (both on the primary image side) for D = 0.5, 0.05, 
and 0.005. We also plot the ratios of (absolute) magni- 
fications of relativistic images of order 1 on the primary 
image side to the secondary image side against M/Dd- 
We finally plot the ratio of magnifications of relativis- 
tic images of orders 1 and 2 (both on the primary im- 
age side) vs M/Dd for same values of V. Compared to 
the ratios of (absolute) magnifications of primary to sec- 
ondary images, the ratios of (absolute) magnifications of 
relativistic images are much less sensitive to changes in 
M/Dd and V. The curves \i\ v /ixi v vs M/Dd, for differ- 
ent values of V, intersect for M/D d « 7.07130 x 10~ 13 . 
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FIG. 9: (color online). Top left: The deflection angles of primary images a p (represented by continuous curves) and of secondary 
images & s (represented by dotted curves) are plotted against M/Dd of MDOs at centers of many galaxies for T> = 0.5,0.05 
and 0.005. Top middle and right: The effective deflection angles of relativistic images of the first order on the same side as 
the primary image df p and on the secondary image side a\ s are plotted against M / Dd of MDOs for the same values of T> as 
in the left figure. Below: The effective deflection angles of relativistic images on the primary image side of the second order 
&2 P and of the first order a\ p vs M/Dd are plotted in the vicinity of zero effective deflection angle for the same values of 
23. The curves for different values of T> on left and right figures intersect, respectively, for (M/Dd)i c ~ 9.33022 x 10~ 13 and 
(M/Dd)ic ~ 9.31854 x 10 -13 . The angular source position (3 — Ifias for all figures. 



For M/Dd less and more than its value on the intersec- 
tion point, (iip/fj,2p are, respectively, higher and lower 
for a lower value of P. As for the case of primary and 
secondary images for any given value of T>, magnifica- 
tions of relativistic images increase with increase in the 
value of M/Dd- However, there is a substantially and 
observationally very important difference in both cases. 
For a given value of M/Dd, just opposite to the case 
of primary and secondary images, the (absolute) mag- 
nifications of relativistic images increase with decrease 
in the value of V. Therefore, sources nearer to the lens 
(with other conditions remaining the same) would give 
relativistic images of higher magnifications. 

In Fig. 11, we first study the variation in the ratio of 
mass of lens to differential time delays among images for 



the change in the value of M/Dd for "D = 0.5, 0.05, and 
0.005. We consider differential time delays between sec- 
ondary and primary images (t s — t p ), the first order rela- 
tivistic image on the primary image side and the primary 
image (ri p — t p ), and second order and first order rela- 
tivistic images both on the primary image side (j2 P — t\ p ). 
We take the angular source position (3 = lfj,as. The ra- 
tio M/(t s — t p ) increases with increase in M/Dd (for a 
fixed value of V) as well as increase in V (for a fixed 
value of M/Dd)- It is obvious from the figure that this 
ratio is very sensitive to changes in distances involved in 
GL. The ratio M/(t\ p — t p ) increases with increase in 
M/Dd (for a fixed value of V); however it decreases with 
increase in T> (for a fixed value of M/Dd). This ratio 
is less (compared to the first we discussed) sensitive to 
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FIG. 10: (color online). Top left and top middle: The magnifications of primary images /i p (represented by continuous curves) 
and the absolute magnifications of secondary images \fi s \ (represented by dotted curves), and their ratios \fj, p /fj, s \ vs M/Dd of 
MDOs at centers of many galaxies are plotted for T> = 0.5, 0.05 and 0.005. Top right and below left: The magnifications of 
relativistic images (on the same side as the primary image) of the first order fii p and of the second order fi2 P are plotted against 
M/Dd of MDOs for same values of T>. Below middle: The ratio of absolute magnifications of relativistic images of first order 
on the primary image side (j,i p to that on the secondary image side \fJ,i B \ are plotted against M/Dd for T> — 0.5, 0.05 and 0.005. 
Below right: The ratio of absolute magnifications of relativistic images (on the primary image side) of the first order fii p to the 
second order /X2 P are plotted against M/Dd for same values of D as in other plots. The curves Hi p /fi2 P vs M/Dd for different 
values of 23 intersect for M/Dd ~ 7.07130 x 10 -13 . For all figures, the angular source position (3 = 1/j.as. 



changes in M/Dd and T). The third ratio M/(t2 P — T\ p ) 
also increases with increase in M/Dd (for a fixed value 
of V); however, for a fixed M/Dd, the variation against 
T> is more fascinating. There exists a critical value of 
M/D d ~ 9.32438 x KT 13 for which curves for different 
values of V intersect. For M/Dd less and more than this 
critical value, the ratio M/(t2 P — Tip) is respectively, less 
and more for smaller value of T>. This ratio is extremely 
insensitive to changes in V and M/Dd, and also in (5 (as 



shown in Section V). Therefore, this ratio can be consid- 
ered almost a constant and can be used to estimate very 
accurate values for masses of lenses once the differential 
time delays T2 P — Ti p are known. We also compute this 
ratio for all MDOs (listed on Table IV) for a few values 
of V > 0.5 (these results are not shown on the plot). 
We find that the slope at any point on M/{j2 P — Tip) 
vs M/Dd curve tends to as V — > 1. We finally plot 
M/(t2 P — Tip) against Dd for different values of V. We 
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FIG. 11: (color online). Top left, top right, and below left: The ratios of mass M of the lens to the differential time delays 
of images vs M/Dd (the ratio of the mass of the lens to its distance) are plotted for T> = 0.5, 0.05 and 0.005. t v and r s 
stand, respectively, for time delays of primary and secondary images, whereas t\ v and T2 P , respectively, represent time delays 
of relativistic images (on the same side as the primary image) of orders 1 and 2. The curves on below left plot intersect for 
M/Dd ~ 9.32438 x 10~ 13 . For a given value of D, ratios of mass to differential time delay are strictly increasing function of 
M/Dd- Below right: The ratios of mass M of the lens to the differential time delay (j2 V — Tip) of images are plotted against 
the distance Dd of the lens for same values of 23. The ratio M/ (j2 V — r\ p ) is not a function of Dd- The angular source position 
j3 = 1 fias for all figures. 



find that the ratio M/(t2 P — T lp ) is extremely insensitive 
to changes in V as well as Dd- Fluctuation in the value of 
this ratio decreases with increase in T>. For a fixed value 
of T>, we find that the ratio M/(t~2 P — Tip) has more than 
1 value for the same value of D^; therefore, this ratio is 
not a function of Dd- However, note that, for any fixed 
value of V, ratios of mass of the lens to the differential 
time delay are strictly increasing functions of M/Dd- 



A. Comparison with Bozza and Mancini's results 

Bozza and Mancini (BM) [l8| obtained differential 
time delays among relativistic images due to GL by a 
general static spherically symmetric spacetime. They 
further modeled MDOs of 12 galaxies as Schwarzschild 
lenses, considered T> — 0.5 (i.e., the lenses to be sym- 
metrically situated between sources and observers) and 
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FIG. 12: (color online). Left: The percentage difference in values for differential time delays (t2 — ri) obtained by Bozza and 
Mancini [1|| (using analytical method) and by us (using numerical method) is plotted against M/Dd of MDOs at centers of a 
few galaxies. t\ and T2, respectively, stand for time delays of relativistic images of orders 1 and 2. Right: The ratios of mass 
M of the lens to the differential time delay (j2 — Ti) vs M/Dd are plotted for both cases. Our numerical results show that the 
ratio M/(t2 — ti) is strictly increasing function of M/Dd- Masses and distances of MDOs are taken from [lj| and are given in 
Table VII in the Appendix. T> — 0.5 and the angular source position j3 = in both figures. 



the angular source position = 0. Using their analytical 
expression for differential time delays among relativistic 
images, they computed differential time delays between 
relativistic images of orders 1 and 2; i.e., (j2 — t\), where 
Ti and T2 are, respectively, time delays of relativistic im- 
ages of orders 1 and 2 for /? = 0. We put their results in 
Table 7 in decreasing order of M/Dd of lenses. We now 
consider the same set of MDOs and use the same values 
for mass M and distance Dd used in their paper. We do 
not use the updated values for M and Dd in this sub- 
section, because we want to compare BM's results with 
ours. Our approach is numerical and we do not take ei- 
ther weak or strong field approximation. As considered 
by those authors, we also take V = 0.5 and (3 = 0, and 
compute differential time delays of relativistic images of 
orders 1 and 2. We then compute percentage difference 
= 100(x — y)/x between ours and their results, where x 
and y are, respectively, differential time delays obtained 
by us and BM. We find that the percentage difference 
ranges approximately from —31.8% to 20.0%, which are 
large. It is possible that MDOs of other galaxies (i.e., 
excluding those considered by BM) give even higher per- 
centage differences. We plot the percentage difference in 
these results against M/Dd (see Fig. 12). Further, using 
BM as well as our results, we compute M/(t2 — Ti) for 
MDOs. We give these results in Table 7. We then plot 
this ratio against M/Dd (see Fig. 12). Results obtained 
by BM show that the ratio M/(t2 — t\) fluctuates quite 
irregularly (showing no rhythm) with an increase in the 
value of M/Dd for otherwise constant situation. There 



seems to be no physical argument in support for this. 
On the other hand, our results show that M/(j2 — Ti) 
is a strictly increasing (though the increase rate is ex- 
tremely small) function of M/Dd- As we did not take 
either weak or strong field approximation at any stage of 
computation and performed numerical computation with 
high precision, our results are very accurate. Therefore, 
we consider percentage differences in results as percent- 
age errors in their results. 

VI. DISCUSSION AND SUMMARY 

It is well-known that the observation of primary and 
secondary images due to GL by an MDO at the center 
of a galaxy is very difficult due to a large extinction of 
electromagnetic radiation (larger extinction for smaller 
wavelength) in the vicinity of a galactic center. In ad- 
dition, radiations at several wavelengths from materials 
accreting on an MDO badly hinders observation of these 
images. These obstacles would be even bigger for rela- 
tivistic images as, compared to primary and secondary 
images, these are formed much closer to the center of 
a galaxy. Unfortunately, observations of relativistic im- 
ages would be much more difficult due to some additional 
reasons. Relativistic images are very much demagnificd, 
unless the lens components (the source, the lens, and the 
observer) are perfectly or highly aligned (/3 << 1/^as), 
and therefore these images are extremely difficult to be 
observed. Supernovae could be more suitable sources 
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for observation of relativistic images, but the probabil- 
ity that a supernova will be highly aligned with the lens 
and observer is extremely small. However, there is a 
silver lining to the demagnification problem associated 
with observation of relativistic images: magnifications of 
relativistic images increase rapidly with the decrease in 
the value of V; i.e., with the decrease in the source-lens 
distance for otherwise constant situation. Thus, sources 
closer to a galactic center would give less demagnified rel- 
ativistic images. Despite this, there is no doubt that the 
observation of relativistic images would be a Herculean 
task. However, with improved observational facilities in 
future and through lucky observations (due to a bright 
source close to a galactic center and highly aligned with 
the galactic center and the observer), the relativistic im- 
ages could possibly be detected some day. Today's these 
almost unthinkable events may be tomorrow's observa- 
tions. The detection of relativistic images would be defi- 
nitely one of the most important discovery in astronomy 
and would have immense implications for general rela- 
tivity and relativistic astrophysics. For examples, these 
observations would provide a test for the general theory 
of relativity in a strong gravitational field. In 0|, we dis- 
cussed that observation of relativistic images would give 
upper bound to the compactness of MDOs and therefore 
would strongly support that these MDOs are black holes. 
The measurements of physical quantities for relativistic 
images would also give very accurate values for masses 
and distances of black holes. 

For any fixed value of V, the ratio of mass M of a 
Schwarzschild lens to differential time delay (t2 P — ti p ) 
or (t2 S — Tig) is not a function of the lens-observer distance 
D<i; however, it is a strictly increasing function of M/Dd. 
Computations show that the ratio M/(j2 P — t\ p ) is ex- 
tremely insensitive to changes in the angular source posi- 
tion as well as the observer- lens and lens-source distances, 
and therefore this awesome physical quantity must be 
treasured as an almost constant for purpose of measure- 
ments. Thus, once we succeed in detecting relativistic 
images and measuring the differential time delay, we can 
immediately compute a very accurate value of mass of the 
MDO acting as a lens. (The accuracy of the result for 
the mass of the MDO will however depend on the accu- 
racy of the measurement of differential time delay.) Our 
computations show that for M/Dd = (0,2.265 x 10~ n ) 
and V= [0.005,1), 

M w 3.734493744773 x 10 5 (r 2p - r lp ), (19) 

where mass M of the MDO and the differential time delay 
i~2 P — T±p are expressed in units of solar mass and minute, 
respectively. As relativistic images can possibly be ob- 
served when the lens components are perfectly or highly 
aligned, we took the angular source position (3 = l/ias 
for computations. However, we found extremely small 
changes in results for computations with /3 = 0. Once, 
the value of mass of the MDO is known, its distance can 
be computed from the results given in Fig. 8 (below). 
Angular separations between relativistic images depend 



on the ratio M/Dd, but fortunately it is extremely insen- 
sitive to the change in the value of V. This would help us 
measure the distance of the MDO very accurately once 
we have the mass of the MDO and the angular separa- 
tion between relativistic images is measured. It is worth 
mentioning that accuracies in determination of distances 
of black holes would, however, depend on accuracies of 
our measurements of differential time delays and angular 
separations between relativistic images. The dependence 
of (absolute) magnifications ratio of relativistic images 
of the first order (i.e., |^i p /^i s |) on M/Dd is extremely 
small (see Fig 10). Therefore, measurement of the mag- 
nifications ratio would give very accurate value for T>. 
This result with already obtained value for observer-lens 
distance Dd would give observer-source distance D s . 

Effective deflection angles of relativistic images play a 
very significant role in analyzing and understanding these 
images. The deflection angles for primary-secondary im- 
age pair as well as relativistic images of Schwarzschild 
black hole lensing are always positive. The effective de- 
flection angles of relativistic images of any order on the 
secondary image side are always positive. However, the 
effective deflection angles of relativistic images of any or- 
der on the primary image side may be positive, zero, or 
negative depending on the value of the angular source 
position (3 and the ratio of mass of the lens to its dis- 
tance (i.e., M/Dd). For a relativistic image (on the pri- 
mary image side) of any order n and for any (nonzero) 
value of M/Dd of the lens, there exists a critical angular 
source position [3 nc such that the effective deflection an- 
gle a^p for that relativistic image is zero. For (3 < (3 nc , 
ot e n p > 0, and for (3 > f3 nc , a e nv < 0. For a given value 
of M/Dd, (3 n c is smaller for smaller n. All sources at 
(3 = (3 nc are lensed to give rise to n th order relativistic 
images (on primary image side) at the same angular po- 
sition 9 np = P„ c . For a fixed value of M/Dd, the angular 
positions of relativistic images are extremely insensitive 
to changes in the angular source position as well as the 
lens-source distance. However, for a theoretical interest, 
it is worth noting that for (3 < (3 nc and (3 > (3 nc , the value 
of 6 np is, respectively, higher and lower for higher value 
of T>. The critical angular source position plays a role 
of flipping point for image positions with respect to the 
change in the value of V. These results help us conclude 
the following: For different sources at the same angular 
position, relativistic images with positive, zero, and neg- 
ative effective deflection angles have, respectively, bigger, 
equal, and smaller (absolute) angular positions for bigger 
values of T>. This is also true for primary and secondary 
images, as deflection angles for them are always positive. 
These conclusions can be also derived from the lens equa- 
tion. Therefore, these results support correctness of our 
numerical computations. 

In GL observations, differential time delays among im- 
ages (not the time delays of individual images) have been 
measured until now. For this reason, studies of time de- 
lays of images have not drawn enough attention. The 
most well-known book on GL gives an expression for time 
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delays of gravitationally lensed images [see Eq. (4.67) in 
||]. The equation has an additive constant term. The 
authors clearly stated that the constant term is the same 
for all rays from the observer to the source plane. Though 
the value for the constant term is not yet determined, the 
expression given in the book is good enough to evaluate 
differential time delays among images. This is because 
the constant term cancels for images of the same source 
by the same lens. Interestingly, Rafikov and Lai [28| have 
recently pointed out that time delays of individual im- 
ages are in fact measurable. This motivated us to first 
compute time delays of images and then use these results 
to compute differential time delays among them. As Eq. 
(4.67) in 0] cannot be used to compute time delays of im- 
ages, we used the method given in Weinberg's book. Our 
results for time delays of primary and secondary images 
turn out to be non-intuitive and very fascinating (see Fig. 
6). For instance, time delays of primary images are al- 
ways smaller for sources nearer to the lens for otherwise 
constant situation. Time delays results for relativistic 
images are also very interesting and important. 

GL as well as gravitational retro-lensing give rise to 
very much demagnified images due to light deflections in 
strong gravitational field. Images due to these two phe- 
nomena can be observationally easily differentiated by 
the fact that the images due to the latter are "orphans" 
in a sense that these are not accompanied by primary- 
secondary images pair as their "parents" . In this paper, 
we have studied only images due to GL. Eiroa and Torres 
[29l | studied retro-lensing by a Schwarzschild black hole. 
They compared magnifications of gravitationally retro- 
lensed images and relativistic images of GL of the same 
order (i.e., the number of turns a light ray makes around 
the lens before reaching the observer) and found that the 
former is significantly greater than the latter. Black holes 
have angular momentum. Therefore, Cunninghom and 
Bardeen (3(| and Rauch and Blandford [3l[ pioneered 
Kerr black hole lensing. As there has been mounting 
observational evidence in support of existence of black 
holes. Kerr lensing has become a very lively research topic 
(see [32[ and references therein) these days. It is worth 
investigating the behavior of the ratio of mass M of the 
lens to differential time delays of images of strong field 
lensing due to changes in (3, M/Dd, T>, and a/m (a is 
the rotational parameter of the Kerr metric). These in- 
vestigations are likely to have immense implication for 
relativistic astrophysics. 

With increasing observational support for MDOs at 
centers of galaxies and stellar size black hole candidates 
to be black holes, the pressure to believe in the existence 
of black holes in the Universe has began to mount. How- 
ever, by the definition of a black hole, there is no and 
there cannot be an iron-clad observational evidence that 
a black hole candidate is indeed a black hole. Given that 
the weak cosmic censorship hyp othesis (WCCH) of Pen- 
rose is still unproven (see [33| and references therein), 
there is no compelling scientific reason to accept that 
all black hole candidates are black holes and none of 



them can be interpreted as a naked (visible) singular- 
ity. Despite the fact that the concept of naked singu- 
larity does not "smell right" to majority of researchers, 
it may not be wise to completely ignore the possibil- 
ity of existence of naked singularities. Researchers think 
that in the vicinity of a spacetime singularity, a myste- 
rious violent marriage of general relativity and quantum 
physics is solemnized and opportunities to observe these 
(through outgoing geodesies from there to us) could help 
us obtaining an unanimously acceptable viable quantum 
gravity theory. Philosophically, it is not clear to us why 
the nature should be malicious to always hide such awe- 
some holy marriages from us. Inspired by these ideas, 
we initiated a new theoretical research project using GL 
phenomena that investigates whether or not black holes 
and naked singularities could be observationally differen- 
tiated (see [lE| Hil HU). Our computations yielded en- 
couraging distinctive results. Whether or not the weak 
cosmic censorship hypothesis of Penrose finally turns out 
to be true, there has to be a cosmic censorship which for- 
bids arbitrary large values of those nakedness parameters 
[e.g., (Q/M) 2 in the Reissner-Nordstrom solution to the 
Einstein-Maxwell equations, where Q and M stand for 
electric charge and mass parameters, respectively] which 
make the system unphysical. Motivated by this idea, 
we hypothesize a new cosmic censorship: Generically, 
marginally and strongly naked singularities do not occur 
in a realistic gravitational collapse. (For definitions of 
weakly naked, marginally naked, and strongly naked sin- 
gularities, see [2l|). The new cosmic censorship hypothe- 
sis (CCH) allows the existence of weakly naked singulari- 
ties, but does not say that these do exist. This hypothesis 
does not imply that the well-known weak cosmic censor- 
ship hypothesis is incorrect. Rather, it says that in case 
the WCCH of Penrose turns out to be incorrect, the new 
cosmic censorship will hold good. In [34| , we showed that 
a Vaidya naked singularity is weakly naked and therefore 
it is not a counter-example to the new CCH. The proof 
of the pudding is in the eating. It may be of an astro- 
physical interest to investigate this subject further. We 
will discuss the new CCH in detail in [35| . 
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VII. APPENDIX 
Tables I- VII 

Tables I through VII are given in this section. 
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TABLE I: Angular positions, bending angles, magnifications, and time delays of primary and secondary images due to GL 
by the Galactic MDO modeled as a Schwarzschild black hole. j3 stands for the angular source position. 9, a, \i, and r 
stand, respectively, for angular positions, deflection angles, magnifications, and time delays of images, with p and s subscripts, 
respectively, for primary and secondary images. t s — t p stands for the differential time delay of the secondary image with respect 
to the primary image. All angles are expressed in arcsec, and time delays and differential time delays are given in minutes. 
(a) The Galactic MDO (lens) has mass M = 3.61 x 1O 6 M , which is at distance D d = 7.62 kpc. M/D d « 2.26 x 10" 11 , where 
M = AlG/c 2 . The ratio of the lens-source distance to the observer-source distance D = 0.5. Results in this table are taken 
from our recent paper [2l| . 
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TABLE II: Effective deflection angles, magnifications, and time delays of relativistic images (on the same side as the primary 
image) due to GL by the Galactic MDO modeled as a Schwarzschild black hole. a e , fj,, and r stand, respectively, for effective 
deflection angles, magnifications, and time delays, with lp and 2p subscripts, respectively, are used for the first and the second 
order relativistic images on the same side as the primary image. T2 P — t\ p , ti p — t p , and t\ p — t s are differential time delays. The 
angular positions of the first and second order relativistic images on the primary image side are, respectively, 8\ p « 24.30283 /ias 
and 02 P ~ 24.27240 fias for all values of angular source position f3 considered in this table, (a) The angular source positions and 
the effective deflection angles are, respectively, expressed in arcsec and fxas, whereas time delays and differential time delays 
are given in minutes. The mass and distance of the lens are as given in (a) of Table I. 
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Mlp 




Tip 




Tip - T s 







48.544793 


X 




48.0274920 


9.666638229 


48.605666 


x 




38.3608537 


23.438763 


23.438763 


10" 


6 


46.544793 


1.34 x 10" 


14 


48.0274920 


9.666638229 


46.605666 


7.21 x 10" 


12 


38.3608537 


23.438764 


23.438762 


10" 


5 


28.544793 


1.34 x 10" 


15 


48.0274920 


9.666638229 


28.605666 


7.21 x 10" 


13 


38.3608537 


23.438771 


23.438754 


10" 


4 


-151.45521 


1.34 x 10" 


16 


48.0274920 


9.666638229 


-151.39433 


7.21 x 10" 


14 


38.3608537 


23.438848 


23.438677 


10" 


3 


-1951.4552 


1.34 x 10" 


17 


48.0274922 


9.666638229 


-1951.3943 


7.21 x 10" 


15 


38.3608540 


23.439616 


23.437910 


10" 


2 


-19951.455 


1.34 x 10" 


18 


48.0275225 


9.666638229 


-19951.394 


7.21 x 10" 


16 


38.3608843 


23.447308 


23.430248 


10" 


1 


-199951.46 


1.34 x 10" 


19 


48.0305628 


9.666638231 


-199951.39 


7.21 x 10" 


17 


38.3639246 


23.525615 


23.354979 


1 




-1999951.5 


1.34 x 10" 


20 


48.3347132 


9.666638248 


-1999951.4 


7.21 x 10" 


18 


38.6680749 


24.463461 


22.721269 



TABLE III: Effective deflection angles, magnifications, and time delays of relativistic images (on the same side as the secondary 
image) due to GL by the Galactic MDO modeled as a Schwarzschild black hole. & e , fi, and r stand, respectively, for effective 
deflection angles, magnifications, and time delays, with subscripts Is and 2s, respectively, are used for the first and the 
second order relativistic images on the secondary image side. ti s — t\ p , t\ s — t p , ti s — r s , and T2 S — t\ p are differential time 
delays. The angular positions of the first and second order relativistic images on the secondary image side are, respectively, 
#i s £s —24.30283 fias and 62s ~ —24.27240 fias for all values of angular source position (3 considered in this table, (a) The 
same as (a) of Table II. 



p 




Mi 




First (outer) relativistic 

Tis Tis " 


ima 
- Tip 


?e 


Tis - T p 


Tis — T s 


"2s 


Second (inner) relativistic ima 

M2s T 2s 


5c 

T2s — Tip 







48.605666 




X 






38.3608537 







23.438763 


23.438763 


48.544793 


X 






48.0274920 


9.666638229 


10" 


6 


50.605666 


-7.21 


x 


10" 


12 


38.3608537 


2.99 x 


10" 


11 


23.438764 


23.438762 


50.544793 


-1.34 x 


10" 


14 


48.0274920 


9.666638229 


10" 


5 


68.605666 


-7.21 


x 


10" 


13 


38.3608537 


2.99 x 


10" 


10 


23.438771 


23.438754 


68.544793 


-1.34 x 


10" 


15 


48.0274920 


9.666638230 


10" 


4 


248.60567 


-7.21 


x 


10" 


14 


38.3608537 


2.99 x 


10" 


9 


23.438848 


23.438677 


248.54479 


-1.34 x 


10" 


16 


48.0274920 


9.666638232 


10" 


3 


2048.6057 


-7.21 


x 


10" 


15 


38.3608540 


2.99 x 


10" 


8 


23.439616 


23.437910 


2048.5448 


-1.34 x 


10" 


17 


48.0274923 


9.666638259 


10" 


2 


20048.606 


-7.21 


x 


10" 


16 


38.3608846 


2.99 x 


10" 


7 


23.447308 


23.430248 


20048.545 


-1.34 x 


10" 


18 


48.0275228 


9.666638528 


10" 


1 


200048.61 


-7.21 


x 


10" 


17 


38.3639276 


2.99 x 


10" 


6 


23.525618 


23.354982 


200048.54 


-1.34 x 


10" 


19 


48.0305658 


9.666641214 


1 




2000048.6 


-7.21 


X 


10" 


18 


38.6681048 


2.99 x 


10" 




24.463491 


22.721299 


2000048.5 


-1.34 x 


10" 


20 


48.3347430 


9.666668077 
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TABLE IV: Masses and distances of MDOs at centers of 40 galaxies are presented in the decreasing order of dimensionless ratio 
of mass to distance [i.e., M/D d = MG[(c 2 D d )] of MDOs. The mass and distance of the Galactic MDO and all other MDOs 
are, respectively, taken from [2(| and [271 ] . 



MDO in 
galaxy 


Mass M 
in Mq 


Distance 
in Mpc 




M 






MDO in 
galaxy 


Mass M 
in M 


Distance 
in Mpc 




M 
D ,i 






Milky Way 


3.61 


X 


10 B 


0.00762 


2 


26467 x 


10" 


11 


NGC5845 


2.4 x 


10 s 


25.9 


4 


42959 x 


10" 


13 


NGC4486(M87) 


3.0 


X 


10 9 


16.1 


8 


90733 x 


10" 


12 


NGC3377 


1.0 x 


10 8 


11.2 


4 


26810 x 


10" 


13 


NGC4649 


2.0 


X 


10 9 


16.8 


5 


69080 x 


10" 


12 


NGC3608 


1.9 x 


10 s 


22.9 


3 


96616 x 


10" 


13 


NGC4594 


1.1 


X 


10° 


9.80 


5 


36561 x 


10" 


12 


NGC4473 


1.1 x 


10 8 


15.7 


3 


34923 x 


10" 


13 


NGC3115 


1.0 


X 


10° 


9.70 


4 


92811 x 


10" 


12 


NGC6251 


5.3 x 


10 8 


93.0 


2 


72424 x 


10" 


13 


NGC224(M31) 


7.0 


X 


10 7 


0.760 


4 


40288 X 


10" 


12 


NGC7052 


3.3 x 


10 8 


58.7 


2 


68737 x 


10" 


13 


IC1459 


2.5 


X 


10 9 


29.2 


4 


09270 x 


10" 


12 


NGC2787 


4.1 x 


10 7 


7.50 


2 


61321 x 


10" 


13 


NGC5128(cenA) 


2.4 


X 


10 s 


4.20 


2 


73158 x 


10" 


12 


NGC4258 


3.9 x 


10 7 


7.20 


2 


58931 x 


10" 


13 


NGC4374(M84) 


1.0 


X 


10 9 


18.4 


2 


59797 x 


10" 


12 


NGC4596 


7.8 x 


10 7 


16.8 


2 


21941 x 


10" 


13 


NGC3998 


5.6 


X 


10 8 


14.1 


1 


89855 x 


10" 


12 


NGC4459 


7.0 x 


10 7 


16.1 


2 


07838 x 


10" 


13 


NGC4486B 


6.0 


X 


10 8 


16.1 


1 


78147 x 


10" 


12 


NGC1023 


4.4 x 


10 7 


11.4 


1 


84502 x 


10" 


13 


NGC4350 


6.0 


X 


10 8 


16.8 


1 


70724 x 


10" 


12 


NGC4564 


5.6 x 


10 7 


15.0 


1 


78463 x 


10" 


13 


NGC4342 


3.1 


X 


10 s 


15.3 


9 


68551 x 


10" 


13 


NGC221(M32) 


2.9 x 


10 6 


0.810 


1 


71145 x 


10" 


13 


NGC3031(M81) 


6.8 


X 


10 7 


3.90 


8 


33483 X 


10" 


13 


NGC821 


8.5 x 


10 7 


24.1 


1 


68599 X 


10" 


13 


NGC4261 


5.2 


X 


10 8 


31.6 


7 


86627 X 


10" 


13 


NGC3384 


1.6 x 


10 7 


11.6 


6 


59348 X 


10" 


14 


NGC4697 


1.7 


X 


10 8 


11.7 


6 


94569 X 


10" 


13 


NGC1068 


1.5 x 


10 7 


15.0 


4 


78027 x 


10" 


14 


CygnusA 


2.9 


X 


10° 


240. 


5 


77616 x 


10" 


13 


NGC4742 


1.4 x 


10 7 


15.5 


4 


31766 x 


10" 


14 


NGC4291 


3.1 


X 


10 s 


26.2 


5 


65604 x 


10" 


13 


NGC7332 


1.5 x 


10 7 


23.0 


3 


11757 x 


10" 


14 


NGC3245 


2.1 


X 


10 s 


20.9 


4 


80314 x 


10" 


13 


NGC2778 


1.4 x 


10 7 


22.9 


2 


92244 x 


10" 


14 


NGC3379 


1.0 


X 


10 8 


10.6 


4 


50969 X 


10" 


13 


NGC4945 


1.4 x 


10° 


3.70 


1 


80875 x 


10" 


14 
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TABLE V: Angular positions, magnifications, and time delays of primary and secondary images due to GL by MDOs (modeled 
as Schwarzschild black holes) at centers of many galaxies. 6, fi, and r stand, respectively, for angular positions, magnifications 
and time delays of images; subscripts p and s attached to them stand, respectively, for primary and secondary images. r s — t v 
stands for the differential time delay of the secondary image with respect to the primary image. The time delays and differential 
time delays are given in minutes, whereas angular positions of images are expressed in arcsec. (a) The first column gives the 
names of galaxies having MDOs with decreasing value of the ratio of mass to the distance (i.e., M/Dd). The ratio of the 
lens-source distance to the observer-source distance T> — 0.5. The angular source position (3 = 1/ias. 



MDO in 

galaxy 







Secondary image 

Ms T s 


r 


3 - Tp 




n 

Up 


Primary 
lip 


image 


Tp 


1\ iT'll TIT 

Milky Way 


— 1 


388176 


-694084.2 


14.92209 


1.71 x 10~ B 


1 


388177 


694085.2 


14.92209 


JN(j-C4486(M87) 


— 


870593 


-435294.8 


12859.74 


Q 


00996 





870594 


435295.8 


12859.74 


t\t r~\ /~i a r a f\ 

NGC4649 


— 


695869 


-347933.4 


8720.112 


(} 


00189 





695870 


347934.4 


8720.110 


NGC4594 


— 


fry r r> nr 

675695 


-337846.2 


4806.677 


o 


00107 





675696 


337847.2 


4806.675 


NGC3115 


— 


647562 


-323779.8 


4383.655 





00101 





647563 


323780.8 


4383.654 


IN Cj0224(M31 ) 


— 


612081 


-306039.8 


308.1496 





00008 





612082 


306040.8 


308. 1495 


IC1459 


— 


590127 


-295062.6 


11035.30 





00278 





590128 


295063.6 


11035.29 


NGC5128(cenA) 


— 


482112 


-241055.1 


1075.303 





00033 





482113 


241056.1 


1075.302 


NGC4374(M84) 


-0 


470173 


-235085.9 


4488.653 





00140 





470174 


235086.9 


4488.652 


NGC3998 


-0 


401930 


-200964.7 


2542.451 





00091 





401931 


200965.7 


2542.450 


NGC4486B 


-0 


389340 


-194669.5 


2730.318 





00101 





389341 


194670.5 


2730.317 


NGC4350 


-0 


381143 


-190570.7 


2734.506 





00103 





381144 


190571.7 


2734.505 


NGC4342 


-0 


287079 


-143539.0 


1441.646 





00071 





287080 


143540.0 


1441.646 


NGC3031(M81) 


-0 


266310 


-133154.8 


317.9070 





00017 





266311 


133155.8 


317.9068 


NGC4261 


-0 


258716 


-129357.8 


2435.988 





00132 





258717 


129358.8 


2435.987 


NGC4697 


-0 


243107 


-121553.1 


799.8507 





00046 





243108 


121554.1 


799.8503 


CygnusA 


-0 


221697 


-110848.0 


13732.21 





00858 





221698 


110849.0 


13732.20 


NGC4291 


-0 


219379 


-109689.4 


1468.994 





00093 





219380 


109690.4 


1468.993 


NGC3245 


-0 


202163 


-101081.3 


1000.754 





00068 





202164 


101082.3 


1000.754 


NGC3379 


-0 


195890 


-97944 


75 


477.5836 





00034 





195891 


97945 


75 


477.5833 


NGC5845 


-0 


194143 


-97071 


07 


1146.906 





00081 





194144 


97072 


07 


1146.905 


NGC3377 


-0 


190571 


-95285 


11 


478.4866 





00034 





190572 


95286 


11 


478.4862 


NGC3608 


-0 


183707 


-91852 


93 


911.4108 





00068 





183708 


91853 


93 


911.4101 


NGC4473 


-0 


168815 


-84407 


35 


530.7089 





00043 





168816 


84408 


35 


530.7085 


NGC6251 


-0 


152252 


-76125 


46 


2575.005 





00228 





152253 


76126 


46 


2575.003 


NGC7052 


-0 


151218 


-75608 


63 


1604.042 





00143 





151219 


75609 


63 


1604.041 


NGC2787 


-0 


149117 


-74558 


08 


199.4783 





00018 





149118 


74559 


08 


199.4781 


NGC4258 


-0 


148433 


-74216 


33 


189.8064 





00017 





148434 


74217 


33 


189.8062 


NGC4596 


-0 


137422 


-68710 


93 


381.5847 





00037 





137423 


68711 


93 


381.5844 


NGC4459 


-0 


132985 


-66491 


96 


343.2016 





00035 





132986 


66492 


96 


343.2012 


NGC1023 


-0 


125296 


-62647 


94 


216.5861 





00023 





125297 


62648 


94 


216.5859 


NGC4564 


-0 


123229 


-61614 


25 


275.9607 





00030 





123230 


61615 


25 


275.9604 


NGC221(M32) 


-0 


120676 


-60337 


77 


14.31073 





00002 





120677 


60338 


77 


14.31072 


NGC821 


-0 


119775 


-59887 


15 


419.6616 





00047 





119776 


59888 


15 


419.6611 


NGC3384 


-0 


074902 


-37450 


84 


81.45874 





00014 





074903 


37451 


84 


81.45860 


NGC1068 


-0 


063777 


-31888 


16 


77.15868 





00015 





063778 


31889 


16 


77.15853 


NGC4742 


-0 


060612 


-30305 


91 


72.24847 





00015 





060613 


30306 


91 


72.24832 


NGC7332 


-0 


051504 


-25751 


92 


78.21023 





00019 





051505 


25752 


92 


78.21004 


NGC2778 


-0 


049866 


-24932 


96 


73.14462 





00018 





049867 


24933 


96 


73.14444 


NGC4945 


-0 


039231 


-19615 


00 


7.424624 





00002 





039232 


19616 


00 


7.424600 



25 



TABLE VI: Angular positions, magnifications, and time delays of the first and second order relativistic images on the same 
side as the primary image. MDOs at centers of many galaxies are modeled as Schwarzschild black hole lenses. 6, fi, and r 
stand, respectively, for angular positions, magnifications, and time delays of images. Subscript p stands for the primary image, 
whereas lp, and 2p stand, respectively, for the first and second order relativistic images on the same side as the primary image. 
Angular positions of the images are expressed in /xas, whereas the time delays and differential time delays are given in minutes. 
Other inputs are the same as (a) of Table V. 



lvi Dyj in 

galaxy 




#2p 


Second order relativistic image 

A*2p 1" 2 p 


T"2p — Tip 






First order relativistic image 
Mlp r lp 


Tip - Tp 


Milky Way 


24.272396 


1 


34 


x 10" 


14 


48.02749 


9.66664 


24.302833 


7 


21 


X 


10" 


12 


38.36085 


23.43876 


NGC4486(M87) 


9 


5467556 


2 


07 


x 10" 


15 


40830.26 


8033.22 


9 


5587269 


1 


12 


X 


10" 


12 


32797.05 


19937.31 


NGC4649 


6 


0993161 


8 


43 


x 10" 


16 


27514.08 


5355.48 


6 


1069644 


4 


55 


X 


10" 


13 


22158.61 


13438.50 


NGC4594 


5 


7507837 


7 


49 


x 10" 


16 


15153.98 


2945.51 


5 


7579950 


4 


05 


x 


10" 


13 


12208.46 


7401.788 


NGC3115 


5 


2818819 


6 


32 


x 10" 


16 


13804.24 


2677.74 


5 


2885053 


3 


41 


x 


10" 


13 


11126.50 


6742.847 


NGC224(M31) 


4 


7189445 


5 


05 


x 10" 


16 


968.8844 


187.442 


4 


7248619 


2 


73 


x 


10" 


13 


781.4426 


473.2931 


IC1459 


4 


3864944 


4 


36 


x 10" 


16 


34662.92 


6694.35 


4 


3919950 


2 


36 


x 


10" 


13 


27968.57 


16933.28 


NGC5128(conA) 


2 


9276717 


1 


94 


x 10" 


16 


3359.469 


642.657 


2 


9313429 


1 


05 


x 


10" 


13 


2716.812 


1641.509 


NGC4374(M84) 


2 


7844704 


1 


76 


x 10" 


16 


14014.24 


2677.74 


2 


7879620 


9 


49 


X 


10" 


14 


11336.50 


6847.846 


NGC3998 


2 


0348357 


9 


38 


x 10" 


17 


7905.583 


1499.53 


2 


0373873 


5 


07 


X 


10" 


14 


6406.049 


3863.599 


NGC4486B 


1 


9093511 


8 


26 


x 10" 


17 


8482.795 


1606.64 


1 


9117454 


4 


46 


X 


10" 


14 


6876.151 


4145.834 


NGC4350 


1 


8297948 


7 


59 


x 10" 


17 


8491.171 


1606.64 


1 


8320893 


4 


10 


X 


10" 


14 


6884.527 


4150.022 


NGC4342 


1 


0380797 


2 


44 


x 10" 


17 


4444.741 


830.099 


1 


0393814 


1 


32 


X 


10" 


14 


3614.642 


2172.996 


NGC3031(M81) 





8933152 


1 


81 


x 10" 


17 


978.3253 


182.086 





8944354 


9 


77 


X 


10" 


15 


796.2390 


478.3322 


NGC4261 





8430953 


1 


61 


x 10" 


17 


7491.180 


1392.42 





8441526 


8 


70 


X 


10" 


15 


6098.755 


3662.769 


NGC4697 





7444293 


1 


26 


x 10" 


17 


2455.980 


455.216 





7453628 


6 


78 


X 


10" 


15 


2000.764 


1200.914 


CygnusA 





6190806 


8 


69 


x 10" 


18 


42071.51 


7765.44 





6198569 


4 


69 


X 


10" 


15 


34306.07 


20573.87 


NGC4291 





6062068 


8 


33 


x 10" 


18 


4499.436 


830.099 





6069670 


4 


50 


X 


10" 


15 


3669.337 


2200.344 


NGC3245 





5147940 


6 


01 


x 10" 


18 


3059.264 


562.325 





5154395 


3 


24 


X 


10" 


15 


2496.939 


1496.185 


NGC3379 





4833420 


5 


2!) 


x 10" 


18 


1458.860 


267.774 





4839481 


2 


86 


X 


10" 


15 


1191.086 


713.5030 


NGC5845 





4747576 


5 


11 


x 10" 


18 


3502.675 


642.657 





4753529 


2 


76 


X 


10" 


15 


2860.018 


1713.113 


NGC3377 





4574487 


4 


74 


x 10" 


18 


1460.666 


267.774 





4580223 


2 


56 


X 


10" 


15 


1192.892 


714.4060 


NGC3608 





4250877 


4 


09 


x 10" 


18 


2779.838 


508.770 





4256208 


2 


21 


X 


10" 


15 


2271.068 


1359.658 


NGC4473 





3589661 


2 


92 


x 10" 


18 


1615.480 


294.551 





3594163 


1 


58 


X 


10" 


15 


1320.929 


790.2203 


NGC6251 





2919802 


1 


93 


x 10" 


18 


7819.582 


1419.20 





2923463 


1 


04 


X 


10" 


15 


6400.381 


3825.378 


NGC7052 





2880290 


1 


88 


x 10" 


18 


4870.271 


883.654 





2883902 


1 


02 


X 


10" 


15 


3986.617 


2382.577 


NGC2787 





2800806 


1 


78 


x 10" 


18 


605.4707 


109.787 





2804318 


9 


60 


X 


10" 


16 


495.6833 


296.2053 


NGC4258 





2775189 


1 


75 


x 10" 


18 


576.0530 


104.432 





2778669 


9 


43 


X 


10" 


16 


471.6212 


281.8150 


NGC4596 





2378733 


1 


28 


x 10" 


18 


1156.050 


208.864 





2381716 


6 


93 


X 


10" 


16 


947.1863 


565.6019 


NGC4459 





2227576 


1 


12 


x 10" 


18 


1038.988 


187.442 





2230370 


6 


07 


X 


10" 


16 


851.5464 


508.3452 


NGC1023 





1977463 


8 


86 


x 10" 


19 


654.7971 


117.821 





1979942 


4 


79 


X 


10" 


16 


536.9766 


320.3907 


NGC4564 





1912746 


8 


2!) 


x 10" 


19 


833.9894 


149.953 





1915144 


4 


48 


X 


10" 


16 


684.0360 


408.0756 


NGC221(M32) 





1834313 


7 


62 


x 10" 


19 


43.22856 


7.76544 





1836613 


4 


12 


X 


10" 


16 


35.46312 


21.15240 


NGC821 





1807017 


7 


40 


x 10" 


19 


1267.462 


227.608 





1809283 


4 


00 


X 


10" 


16 


1039.854 


620.1931 


NGC3384 





0706679 


1 


13 


x 10" 


19 


243.5083 


42.8438 





0707566 


6 


11 


X 


10" 


17 


200.6645 


119.2059 


NGC1068 





0512343 


5 


95 


x 10" 


20 


229.8712 


40.1661 





0512985 


3 


21 


X 


10" 


17 


189.7051 


112.5466 


NGC4742 





0462761 


4 


85 


x 10" 


20 


215.0139 


37.4883 





0463341 


2 


62 


X 


10" 


17 


177.5255 


105.2772 


NGC7332 





0334136 


2 


53 


x 10" 


20 


231.9743 


40.1661 





0334555 


1 


37 


X 


10" 


17 


191.8082 


113.5982 


NGC2778 





0313223 


2 


22 


x 10" 


20 


216.8061 


37.4883 





0313615 


1 


20 


X 


10" 


17 


179.3178 


106.1734 


NGC4945 





0193859 


8 


52 


x 10" 


21 


21.90093 


3.74883 





0194102 


4 


60 


X 


10" 


18 


18.15210 


10.72750 



TABLE VII: Comparison of Bozza and Mancini's approximate analytical pl| and our numerical results. Mass and distance of 
MDOs are given in the units of solar mass and Mpc, respectively. M/Dd in the third column is dimensionless (M = MG/c 2 ). 
Differential time delays T2 — ri are expressed in minutes. T2 and n stand, respectively, for time delays of the first and second 
order relativistic images for the angular source position j3 = 0. The ratio of source-lens to source-observer distances T> = 0.5. 
The ratio of mass of an MDO to the differential time delay M/(t2 — ri) is expressed in terms of solar mass/minute. Masses 
and distances of MDOs in this table are taken the same as in I Lai. 



MDO in 

galaxy 


Mass 
M 


Distance 
D d 


M/D d 


Differential time delay (t2 — ti) 
Analytical Numerical % difference 


M/(T2 

Analytical 


-ti) 


Numerical 


Milky Way 


2 


8 X 10° 


0.0085 


1 


57 > 


; 10" 


ii 


6 


7.4977 


20.0 


466666.666666667 


373449 


374476538 


NGC3115 


2 


x 10 9 


8.4 


1 


14 > 


; 10" 


n 


5430 


5355.5 


-1.4 


368324.125230203 


373449 


374476535 


NGC4486(M87) 


3 


3 x 10 9 


15.3 


1 


03 > 


; 10" 


n 


8958 


8836.5 


-1.4 


368385.800401875 


373449 


374476534 


NGC4594 


1 


x 10 9 


9.2 


5 


20 > 


< 10" 


12 


2712 


2677.7 


-1.3 


368731.563421829 


373449 


374476530 


NGC4374 (M84) 


1 


4 x 10 9 


15.3 


4 


37 > 


< 10" 


12 


3798 


3748.8 


-1.3 


368615.060558189 


373449 


374476529 


NGC224(M31) 


3 


x 10 7 


0.7 


2 


05 > 


< 10" 


12 


84 


80.332 


-4.6 


357142.857142857 


373449 


374476527 


NGC4486B(M104) 


5 


7 x 10 s 


15.3 


1 


78 > 


< 10" 


12 


1548 


1526.3 


-1.4 


368217.054263566 


373449 


374476527 


NGC4342(IC3256) 


3 


x 10 s 


15.3 


9 


37 > 


; 10" 


13 


816 


803.32 


-1.6 


367647.058823529 


373449 


374476527 


NGC3377 


1 


8 x 10 s 


9.9 


8 


69 > 


; 10" 


13 


486 


481.99 


-0.8 


370370.370370370 


373449 


374476527 


NGC4261 


4 


5 x 10 s 


27.4 


7 


85 > 


; 10" 


13 


1224 


1205.0 


-1.6 


367647.058823529 


373449 


374476526 


NGC7052 


3 


3 x 10 s 


58.7 


2 


69 > 


; 10" 


13 


894 


883.65 


-1.2 


369127.516778523 


373449 


374476526 


NGC0221(M32) 


3 


4 x 10 6 


0.7 


2 


32 > 


; 10" 


13 


12 


9.1043 


-31.8 


283333.333333333 


373449 


374476526 



20 



[1] In gravitational lensing theory, weak and strong fields 
usually refer to the regions of spacetime which cause, 
respectively, small and large deflections of light. 

[2] A. Eddington, Space, Time and Gravitation (Cambridge 
Univ. Press, Cambridge, 1920). 

[3] P. Schneider, J. Ehlers, and E. E. Falco, Gravitational 
Lenses ( Springer- Ver lag, Berlin, 1992). 

[4] P. Schneider, C. Kochanek, and J. Wambsganss, Grav- 
itational Lensing: Strong, Weak and Micro (Springer- 
Verlag, Berlin, 2006). 

[5] R. Narayan and M. Bartelmann, |astro-p h/9606001 
O . Wucknitz , i arXiv:0709.4005l [astro-ph]^ 

[6] C. Darwin, Proc. R. Soc. London A 249, 180 (1959); 
263, 39 (1961). 

[7] R. Atkinson, The Astronomical Journal, 70, 517 (1965). 

[8] S. Frittelli and E. T. Newman, Phys. Rev. D 59, 124001 
(1999). 

[9] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D 62, 
084003 (2000). 

[10] S. Frittelli, T. P. Kling, and E. T. Newman, Phys. Rev. 
D 61, 064021 (2000); T. P. Kling, E. T. Newman, and 
A. Perez, ibid 62, 024025 (2000);62, 109901(E) (2000). 

[11] V. Perlick, larXiv:0708.0178l [gr-qc] . 

[12] V. Perlick, Living Rev. Relativity 7, 9 (2004); Phys. 
Rev. D 69, 064017 (2004); Commun. Math. Phys. 220, 
403 (2001); W. Hasse and V. Perlick, J. Math. Phys. 
47, 042503 (2006). W. Hasse and V. Perlick, Gen. Rela- 
tiv. Gravit. 34, 415 (2002); T. Foertsch, W. Hasse, and 
V. Perlick, Classical Quantum Gravity 20, 4635 (2003). 

[13] P. Amore, M. Cervantes, A. De Pace, and F. M. Fer- 
nandez, Phys. Rev. D 75, 083005 (2007); P. Amore and 
S. A. Diaz, ibid 74, 083004 (2006); P. Amore, S. Arceo, 
and F. M. Fernandez, ibid 73, 083004 (2006). 

[14] A. F. Zakharov, F. De Paolis, G. Ingrosso, and 
A. A. Nucita, Astron. Astrophys. 442, 795 (2005); New 
Astron. 10, 479 (2005); A. F. Zakharov and Yu. V. Bary- 
shev, Int. J. Mod. Phys. D 11, 1067 (2002); Classical 
Quantum Gravity 19, 1361 (2002). 

[15] G. N. Gyulchev and S. S. Yazadjiev, Phys. Rev. D 
75, 023006 (2007); R. Whisker, ibid 71, 064004 (2005); 

C. R. Keeton and A. O. Potters, ibid 73, 044024 (2006); 
73, 104032 (2006); E. F. Eiroa, G. E. Romero, and 

D. F. Torres, ibid 66, 024010 (2002); E. F. Eiroa, ibid 
73, 043002 (2006); A. O. Petters, Mon. Not. R. As- 
tron. Soc. 338, 457 (2003); S. V. Iyer and A. O. Pet- 
ters, Gen. Relativ. Gravit. 39, 1563 (2007); N. Mukher- 
jee and A. S. Majumdar, ibid 39, 583 (2007); A. S. Ma- 
jumdar and N. Mukherjee, Mod. Phys. Lett. A 20, 2487 
(2005); E. F. Eiroa, Braz. J. Phys. 35, 1113 (2005); 

E. F. Eiroa, G. E. Romero, Phys. Letts. B663, 377 
(2008); G. S. Bisnovatyi-Kogan and O. Yu. Tsupkol, 
Astrophysics, 51 99 (2008). 

[16] S. Fernando and S. Roberts, Gen. Relativ. Gravit. 34, 

1221 (2002); P. Mutka and P. Mahonen, Astrophys. J. 

581 1328 (2002); C. Stornaiolo, Gen. Relativ. Gravit. 

34, 2089 (2002); P. Mutka and P. Mahonen, Astrophys. 

J. 581 1328 (2002); K. Lake, Phys. Rev. D 65, 087301 

(2002); A. Bhadra, Phys. Rev. D 67, 103009 (2003); 

R. A. Konoply a, Phys. Letts. B 644, 219 (2007); J. Briet 

and D. Hobill. larXiv:0801. 3859 [gr-qc]. 
[17] V. Bozza, S. Capozziello, G. Iovane, and G. Scarpetta, 



Gen. Relativ. Gravit. 33, 1535 (2001). 
[18] V. Bozza and L. Mancini, Gen. Relativ. Gravit. 36, 435 
(2004). 

[19] K. S. Virbhadra, D. Narasimha, and S. M. Chitre, As- 
tron. Astrophys. 337, 1 (1998). 

[20] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D 65, 
103004 (2002). 

[21] K. S. Virbhadra and C. R. Keeton, Phys . Rev. D 77, 

124014 (2008). See also larXiv:0710.2333l V r l [gr-qc] for 

more results in tables. 
[22] G. N. Gyulchev, S. S. Yazadjiev. larXiv:0806.3289l [gr-qc] ; 

M. C. Werner and A. O. Petters, Phys. Rev. D 76, 064024 

(2007). 

[23] M. Safonova and D. F. Torres, Mod. Phys. Lett. A 17, 
1685 (2002); J. M. Tejeiro S. and E. A. Larranaga R., 
|gr-qc/0505054} K. K. Nandi, Phys. Rev. D 74, 024020 
(2006); F. Rahaman, M. Kalam, and S. Chakraborty, 
Chin. J. Phys. 45, 518 (2007); T. K. Dey and S. Sen, 
Mod. Phys. Lett. A A23, 953 (2008). 

[24] V. Kagramanova, J. Kunz, and C. Lammerzahl, Phys. 
Letts. B634, 465 (2006); I. Ciufolini and F. Ricci, Clas- 
sical Quantum Gravity 19, 3863 (2002); 19, 3875 (2002); 
S. S. Yazadjiev, |gr-qc/0 305053 Y. Chen, Z. Chen, 
and J. Wang, Gen. Relativ. Gravit. 39, 401 (2007); 
M. P. Dabrowski and F. E. Schunck, Astrophys. J. 535, 
316 (2000); T. Matos and R. Becerril, Classical Quantum 
Gravity 18, 2015 (2001); V. A. De Lorenci, N. Figueiredo, 
H. H. Fliche, and M. Novello, Astron. Astrophys. 369, 
690 (2001); F. E. Schunck, B. Fuchs, and E. W. Mielke, 
Mon. Not. R. Astron. Soc. 369, 485 (2006); K. Sarkar 
and A. Bhadra, Classical Quantum Gravity 23, 6101 
(2006); Y. Chen, Z. Chen, and J. Wang, Gen. Relativ. 
Gravit. 39, 401 (2007). 

[25] S. Weinberg, Gravitation and Cosmology: Principles and 
Applications of the General Theory of Relativity (Wiley, 
New York, 1972). 

[26] F. Eisenhauer et al, Astrophys. J. 6 28, 246 (2005). 

[27] K. Gebhardt, |http://hoku.as.utexas.edu/^gebhardt 
/blackhole.html. 

[28] R. R. Rafikov and D. Lai, Phys. Rev. D 73, 063003 
(2006). 

[29] E. F. Eiroa and D. F. Torres, Phys. Rev. D 69, 063004 
(2004). 

[30] C. T. Cunningham and J. M. Bardeen, Astrophys. J., 

183, 237 (1973). 
[31] K. P. Rauch and R. D. Blandford, Astrophys. J. 421, 46 

(1994). 

[32] S. E. Vazquez and E. P. Esteban, Nuovo Cim. B 119, 
489 (2004). 

[33] R. Penrose, in Black Holes and Relativistic Stars, edited 
by R. M. Wald (The University of Chicago Press, 
Chicago, 1998), p. 103; K. S. Virbhadra, S. Jhingan, and 
P. S. Joshi, Int. J. Mod. Phys. D 6, 357 (1997). K. S. Virb- 
hadra, Int. J. Mod. Phys. A 12, 4831 (1997); Phys. Rev. 
D 60, 104041 (1999). 

[34] C-M. Claudel, K. S. Virbhadra, and G. F. R. Ellis, J. 
Math. Phys. 42, 818 (2001). 

[35] K. S. Virbhadra, Gravitational lensing by Reissner- 
Nordstrom black holes and naked singularities (in prepa- 
ration) . 



